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ABSTRACT 


A  three-dimensional  objective  analysis  technique  known  as  CRAM  (Conditional 
Relaxation  Analysis  Method)  was  applied  to  investigate  various  properties  of  ballistic 
winds  on  a  mesoscale  in  mountainous  regions.  From  a  12-day  sample  of  upper-air 
soundings  taken  5  times  a  day  at  2-hr  intervals  for  12  rawinsonde  stations  in  the 
Ft.  Huachuca  region  of  southeastern  Arizona,  and  artillery  firings  taken  twice  a  day, 
CRAM  analyses  oi  temperature,  density,  and  winds  were  performed  for  10  atmospheric 
zones  between  the  surface  and  8,000  m  using  an  IBM-7094. 

It  was  determined  that  the  CRAM  technique  produced  fields  which  fpd  the  desir¬ 
able  features  of  map  winds,  i.e.,  the  contour  patterns  were  relatively  smiopth  and  varied 
slowly  with  time.  The  residual  deflection  errors  which  resulted  were  smaller  for 
CRAM  (75.2  m)  than  for  a  single  station  (Ft.  Huachuca)  near  the  firing  range  (85.1  m). 

It  was  also  found  that  the  time  decay  of  ballistic  winds  in  the  firing  area  was  smaller 
using  CRAM  than  using  the  Ft.  Huachuca  observation,  which  implies  that  CRAM  is  a 
better  tool  with  which  to  make  a  persistence  forecast  than  a  single  station. 
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1.0  INTRODUCTION 


In  the  general  artillery  firing  problem  it  is  necessary  for  the  gunner  to  apply 
various  corrections  so  that  a  round,  or  series  of  rounds,  from  his  piece  will  come  as 
near  as  possible  to  striking  a  preassigned  target.  There  are  two  major  factors  in 
artillery  correction  problems.  The  first  has  to  do  mainly  with  characteristics  of  the 
weapon  system  being  employed.  These  include  projectile  weight,  powder  temperature, 
muzzle  velocity,  and  drift,  and  their  investigation  is  essentially  beyond  the  scope  of 
this  project.  The  second  factor  to  be  considered  is  the  influence  of  the  atmosphere 
through  which  the  projectile  travels.  Of  Importance  here  are  wind,  temperature,  and 
density.  Although  the  method  of  registration  (essentially  a  trial  and  error  method  of 
correcting  subsequent  shots  based  on  results  of  previous  ones,  while  firing  a  series 
of  rounds)  is  possible  in  principle,  it  is  not  the  most  desirable  tactically  (the  element 
of  surprise  may  be  of  importance),  practically  (registration  requires  visual  observa¬ 
tion  of  the  impact),  or  economically  (it  can  be  costly). 

At  present,  the  U.S.  Army  makes  ballistic  corrections  for  meteorological  effects 
on  an  artillery  projectile  based  on  a  single  sounding  method.  An  Artillery  Metro 
Section  launches  a  radiosonde  and  from  the  data  computes  a  ballistic  message  which  is 
provided  to  the  gunner.  This  observation  may  be  removed  in  both  time  and  space  from 
the  projectile  firing.  Research  results  \Lowenthal,  1953,  1957)  indicate  that  improved 
ballistic  corrections  can  be  obtained  by  use  of  multi-sounding  data  from  which  fields 
of  wind,  temperature,  and  density  can  be  mapped.  This  so-called  “map  wind  technique” 
has  the  desirable  property  that  the  analyzed  fields  arp  fairly  smooth  and  slowly  varying 
in  time,  and  as  such  present  a  hopeful  avenue  for  prognosis. 

The  objective  of  this  program  was  to  develop  an  automated  objective  analysis 
technique  with  emphasis  on  the  mesoscale  meteorological  features  to  produce  map 
winds,  temperatures,  and  densities,  which,  when  incorporated  in  an  integrated  metro 
message,  would  minimize  the  error  from  artillery  firings  made  at  the  same  time. 

The  research  was  to  be  conducted  with  a  view  toward  providing  the  best  possible 
answer  to  the  following  questions: 

(a)  To  what  extent  can  artillery  corrections  be  improved  by  using 
integrated  map  winds  instead  of  a  single  sounding? 
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(b)  How  does  mountainous  terrain  affect  the  variability  of  ballistic 
winds? 

(c)  What  is  the  optimum  placement  and  spacing  of  artillery  metro 
sections  in  mountainous  terrain? 

(d)  How  often  should  soundings  be  made  in  mountainous  terrain? 

(e)  How  accurately,  and  for  what  duration,  can  map  winds  be 
forecast? 

For  this  program,  a  field  experiment  was  conducted  at  Ft.  Huachuca  in  which 
multiple  sounding  data  were  taken  concurrently  with  artillery  firings  for  a  number  of 
days  in  January  and  February  1965.  A  detailed  description  of  the  field  experiment  is 
given  in  Section  2.0. 

An  analysis  technique  suited  to  the  ballistic  wind  problem  (called  CRAM— 
Conditional  Relaxation  Analysis  Method)  was  chosen  as  the  primary  investigatory  tool. 
It  is  a  technique  that  the  contractor  has  had  previous  experience  with  and  was  extend¬ 
able  to  three  dimensions  with  relatively  minor  modifications.  The  technique  is 
described  in  Section  3.0,  The  detailed  program  descriptions  are  contained  in  the 
supplement  to  this  report. 

While  the  advantages  of  map  ballistic  winds  over  single  station  ballistic  winds, 
in  principle ,  are  readily  apparent  (cf.  Lowenthal,  1953),  an  important  evaluation 
criterion  is  how  well  a  given  technique  performs  when  used  for  correcting  actual 
artillery  firings.  This  is  discussed  in  Sections  4.0  and  7.0. 

Questions  concerning  characteristics  of  the  winds  themselves,  such  as  time  and 
space  variability  and  predictability,  are  discussed  in  Sections  5.0,  6.0,  and  8.0. 
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2.0  DESIGN  OF  THE  FIELD  EXPERIMENT 


Artillery  firings  were  conducted  at  Ft.  Huachuca,  Arizona  on  12  days  during 
January  and  February^  1965  as  part  of  this  study.  Two  8-inch  howitzers  were  fired 
at  1000  MST  and  1400  MST  at  targets  approximately  14,500  m  east-northeast  from  the 
gun  locations.  Each  piece  fired  ten  rounds  alternately  at  one-minute  intervals  using 
fixed  quadrant  elevation  angles  and  azimuths  (see  Table  2-1).  Impact  points  were 
determined  visually  and  the  center-of-impact  (Cl)  for  the  ten  rounds  from  each  piece 
was  computed.  The  maximum  ordinate  of  the  projectiles  was  approximately  9700  m, 
or  about  8200  m  above  the  gun  level. 

Meteorology's!  data  were  collected  from  a  network  of  twelve  rawinsonde  stations 
surrounding  the  firing  range  (see  Fig.  2-1).  Soundings  were  taken  5  times  a  day  at  two 
hour  intervals  beginning  at  0600  MST,  for  even  numbered  days  January  2—12,  odd 
numbered  days  January  25—31,  and  even  numbered  days  February  2—12,  1965.  Thus, 
there  were  16  days  on  which  soundings  were  taken,  or  80  observation  times  in  all. 

The  meteorological  data  were  obtained  using  AN/GMD-1A  equipment.  The  twelve 
meteorological  stations  for  which  rawinsonde  data  were  collected  are  listed  in  Table 
2-2. 

The  meteorological  data  were  punched  on  cards  and  transferred  to  magnetic 
tape  by  the  sponsoring  agency  using  a  Burroughs  B-5000  Computer. 


TABLE  2-1 
FIRING  DATA 


Gun  no. 

Location* 

Quadrant  elev. 
(mils) 

Easting 

Northing 

Azimuth  (mils) 

1 

555472.8 

3493987.7 

1325.2 

1114 

2 

555469.4 

3494002.5 

1322.8 

1117 

♦Coordinates  are  based  on  the  Universal  Transverse  Mercator 
Grid  (Zone  12). 
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TABLE  2-2 

METEOROLOGICAL  STATIONS  IN  THE  FT.  HUACHUCA  NETWORK 


Station 

Easting(I)* 

Northing(J)* 

Elevation(m) 

Benson 

57.01 

353.59 

1083 

Bisbee 

60.68 

347.01 

1445 

Douglas 

64.10 

346.82 

1255 

Fairbanks 

56.81 

350.92 

1282 

Hereford 

58,62 

347.72 

1290 

Ft.  Huachuca 

56.25 

349.36 

1433 

Nogales 

51.48 

347.51 

1190 

Patagonia 

52.40 

348.95 

1240 

Parker  Canyon 

55.37 

347.51 

1685 

Sonoita 

53.29 

350.41 

1480 

Tombstone 

59.26 

350.38 

1450 

Wilcox 

60.49 

356.73 

1275 

♦Grid  coordinates  are  based  on  the  10,000-meter  Universal 
Transverse  Mercator  Grid  (Zone  12). 


3.0  DATA  PROCESSING 


The  meteorological  data  described  in  the  previous  section  needed  considerable 
editing,  error  checking,  calculation  of  additional  variables,  etc.,  before  objective 
analysis  methods  could  be  employed.  The  procedures  used  were  designed  to  be  as 
consistent  as  possible  with  U.S.  Army  procedures  for  performing  ballistic  computa¬ 
tions.  The  major  difference  between  current  practices  and  those  adopted  in  this  study 
was  the  referencing  of  meteorological  zone  data  to  a  geographical  location  representa¬ 
tive  of  the  radiosonde’s  actual  location,  rather  than  the  customary  practice  of  assigning 
the  launching  station’s  geographical  coordinates  to  a  piece  of  meteorological  informa¬ 
tion  regardless  of  how  far  away  from  that  station  the  balloon  had  subsequently  drifted. 
Also  taken  into  account  in  this  study  was  the  referencing  of  the  balloon’s  height  to  the 
terrain  over  which  it  was  passing,  rather  than  to  the  elevation  of  the  launch  station. 

3.1  Computation  of  Ballistic  Quantities 

3.1.1  Conventional  Methods 

Ballistic  wind,  temperature,  and  density  are  computed  in  order  to  apply  artillery 
corrections  using  known  unit  effects  from  firing  tables.  To  arrive  at  these  ballistic 
quantities  for  a  particular  station,  measurements  of  upper-air  parameters  are  made 
by  means  of  a  radiosonde  which  is  tracked  by  a  rawin  set  (AN/GMD-1).  The  location 
of  the  balloon  at  each  artillery  zone  limit  (see  Table  3-1),  as  projected  to  the  earth’s 
surface,  is  plotted  on  a  plotting  board.  From  these  plots,  the  average  wind  speed  and 
direction  for  each  of  the  atmospheric  zones  is  determined.  The  computation  of  these 
zone  winds  is  a  preliminary  step  in  the  determination  of  ballistic  winds.  The  mean 
virtual  temperature  and  mean  density  for  each  zone  are  also  determined  and  pre¬ 
sented  in  terms  of  percent  of  standard  as  referenced  to  mean  standard  zone  density 
and  temperature  based  on  the  ICAO  (International  Civil  Aviation  Organization)  atmos¬ 
phere.  The  zone  values  of  wind,  density,  and  temperature  are  then  weighted  according 
to  specified  zone  weighting  factors,  with  the  resulting  weighted  quantities  being  the 
ballistic  values.  The  weighting  factors,  developed  by  ballisticians,  are  used  to 
establish  the  proportional  effect  of  the  meteorological  conditions  in  each  zone  upon 
the  total  effect  exerted  by  the  atmosphere  through  which  a  projectile  passes.  Table 
3-1  shows  the  weighting  factors  for  the  artillery  problem  in  the  Ft.  Huachuca  experi¬ 
ment  where  the  maximum  ordinate  of  the  trajectory  was  in  Zone  10  (these  are  known 
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TABLE  3-1 

ZONE  STRUCTURE  AND  WEIGHTING  FACTORS  FOR  BALLISTIC 
COMPUTATIONS  (maximum  ordinate:  8,000  m) 


Zone 

Height  (m) 

Weighting  factors* 

Base 

Top 

Wind 

Temperature 

Density 

1 

Surface 

200 

0.01 

0.00 

0.03 

2 

200 

500 

0.02 

0.01 

0.04 

3 

500 

1000 

0.02 

0.02 

0.07 

4 

1000 

1500 

0.04 

0.02 

0.07 

5 

1500 

2000 

0.03 

0.03 

0.07 

6 

2000 

3000 

0.07 

0.07 

0.13 

7 

3000 

4000 

0.08 

0.13 

0.12 

8 

4000 

5000 

0.09 

0.22 

0.11 

9 

5000 

6000 

0.09 

0.23 

0.11 

10 

6000 

8000 

0.55 

0.27 

0.25 

♦From  FM  6-16,  Tables  for  Artillery  Meteorology. 


as  “Line  10”  weighting  factors).  The  zone  values  for  Zones  1,  2,  3,  ...,  10  are 
weighted  and  summed  to  arrive  at  the  ballistic  values  for  Line  10.  When  there  is  a 
difference  in  altitude  between  the  artillery  piece  and  the  meteorological  station,  the 
temperature  and  density  values  are  corrected  using  calculations  based  on  the  ICAO 
standard  atmosphere.  No  effort  is  made  to  adjust  the  ballistic  winds  for  such  a  height 
difference  because  there  is  no  known  specific  relation  between  the  speed  and  direction 
of  the  wind  and  this  height  difference. 

3.1.2  Modified  Methods 

Objective  analyses  of  individual  atmospheric  zones  are  required  to  construct 
analyzed  fields  of  ballistic  wind,  temperature,  and  density.  In  preparing  the  data  for 
analysis,  it  was  decided  to  account  for  balloon  drift  such  that  data  values  for  individual 
zones  would  be  assigned  horizontal  coordinates  based  on  the  location  of  the  balloon  for 
the  midpoint  of  the  zone.  When  dealing  with  conventional  synoptic-scale  data  it  is 
customary  to  ignore  balloon  drift  with  altitude.  However,  for  a  mesoscale  problem, 
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as  pointed  exit  by  Fujita  (1960),  a  significant  error  could  result  in  the  analyses  of 
higher  zones  if  it  were  assumed  that  the  observations  were  representative  of  points 
directly  above  the  launch  sites.  Figure  3-1  shows  the  displaced  locations  of  the 
several  radiosondes  from  their  launch  sites  for  the  Ft.  Huachuca  network  at  1400  MST 
on  January  25,  1965  upon  reaching  the  mid-level  of  Zone  10.  Of  particular  importance 
is  the  fact  that  Zone  10  receives  the  largest  weight  in  computing  Line  10  ballistic 
winds.  Note  that  all  the  balloons  have  moved  to  the  east  or  southeast  of  Ft.  Huachuca. 
The  radiosonde  launched  at  Ft.  Huachuca  is  actually  located  nearly  10  km  southeast  of 
Hereford  at  Zone  10. 

Another  consideration  in  analyzing  data  from  a  multi-sounding  network  has  to  do 
with  the  manner  in  which  height-above-terrain  is  handled.  As  stated  above,  zone  winds 
are  computed  using  the  height  above  terrain  rather  than  referencing  observations  to 
the  height  of  the  artillery  piece.  In  this  study,  the  terrain  height  at  the  location  of  the 
balloon,  and  not  the  terrain  height  at  the  launch  site,  was  used  in  determining  zone 
heights.  This  distinction  can  be  significant  in  mountainous  terrain,  particularly  for 
the  lower  zones. 

3.2  Preprocessing  the  Data 

The  meteorological  data  contained  on  magnetic  tapes  represent  approximately 
900  radiosonde  flights  consisting  of  roughly  20  significant  data  levels  each  and  50  or 
so  detailed  wind  data  levels  each,  or  about  240,000  relevant  numbers  on  tape.  With 
this  quantity  of  data  involved,  certain  gross  error  checking  procedures  are  recom¬ 
mended.  With  the  aid  of  the  computer,  the  data  were  checked  for  such  things  as 
proper  headings,  station  elevations,  dates  and  times,  etc.  Routines  which  checked  to 
assure  that  times  and  heights  increased  and  pressures  decreased  monotonically  were 
useful  in  uncovering  bad  data  or  data  out  of  sequence.  An  additional  test  was  made  of 
times  and  heights  based  on  the  expected  rate  of  balloon  ascent.  When  this  phase  was 
completed,  a  relatively  error-free  meteorological  data  sample  was  available  for 
additional  processing. 

The  next  phase  consisted  of  the  computation  of  zone  winds,  temperatures,  and 
density  from  the  raw  data  for  each  radiosonde  flight.  The  zone  winds  were  computed 
using  the  time  and  horizontal  location  at  which  the  balloon  entered  the  base  of  a 
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particular  zone  and  similar  information  when  the  balloon  reached  the  top  of  the  zone. 
The  mean  virtual  temperature  for  a  zone  was  calculated  from  the  significant  level 
temperature  and  humidity  data.  Density  was  calculated  from  the  equation  of  state 
using  the  mean  virtual  zone  temperature  and  the  pressure  at  the  midpoint  of  the  zone. 
Temperatures  and  densities  were  referenced  to  the  ICAO  standard  atmosphere  and 
expressed  in  terms  of  departure  from  standard  in  percent.  The  data  were  then  ready 
for  objective  analysis. 

3.3  The  CRAM  Objective  Analysis  Procedure 

One  objective  of  this  study  is  to  develop  an  automated  objective  analysis  tech¬ 
nique  to  produce  map  wind  fields  from  which  ballistic  quantities  may  be  estimated. 

Objective  analysis  is  concerned  with  defining  the  spatial  distribution  of  a  vari¬ 
able  at  a  set  of  regularly  spaced  grid  points  given  an  irregularly  spaced  set  of 
observations.  The  technique  chosen  for  this  study  is  one  based  on  a  generalized  version 
of  Carstensen’s  relaxation  method  (1962),  which  was  developed  and  applied  by  Thomasell 
and  Welsh  (1963),  called  CRAM  (the  Conditional  Relaxation  Analysis  Method).  It  is 
ideally  suited  to  a  three-dimensional  analysis  problem  such  as  this  one  and  is  relatively 
fast  from  a  computer  time  consideration. 

In  CRAM,  the  procedure  for  interpolating  between  observations  requires  that  the 
analyzed  grid-point  values  satisfy  Poisson’s  equation,  subject  to  the  constraints  imposed 
by  tne  observations  and  an  arbitrarily  defined  set  of  boundary  values.  The  observations 
determine  the  analysis  through  their  role  as  internal  boundary  points  in  the  solution  of 
Poisson’s  equation.  Poisson’s  equation  in  finite -difference  form  is  solved  numerically 
by  a  relaxation  procedure. 

3.3.1  The  Analysis  Grid 

The  analysis  grid  is  based  on  the  Universal  Transverse  Mercator  Grid  system 
and  is  shown  in  Fig.  2-1.  The  horizontal  grid  spacing  is  10  km.  The  I  and  J  directions 
of  the  analysis  grid  are  essentially  east-west  and  north-south  respectively.  The 
analysis  area  consists  of  17  grid  units  in  the  I-direction  (from  I  =  50  to  I  =  66)  and 
15  grid  units  in  the  J-direction  (from  J  =  345  to  J  =  359)  or  about  170  km  x  150  km. 

The  vertical  coordinate  (K)  consists  of  the  10  atmospheric  zones  and  is  of  variable 
spacing,  being  250  m  between  Zones  1  and  2  and  1500  m  between  Zones  9  and  10. 
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Two  additional  zones — Zone  “0"  and  Zone  11— are  added  to  supply  boundary  values  at 
the  edges  of  the  analysis  volume. 

Subsequent  sections  describe  the  various  major  components  of  CRAM:  the  initial 
guess  technique,  correction  procedures,  relaxation  method,  smoothing,  and  ballistic 
computation.  Detailed  program  descriptions  are  contained  in  the  supplement  to  this 
report. 

3.3.2  _  The  Initial  Guess 

CRAM  requires  the  use  of  an  initial-guess  field,  which  undergoes  subsequent 
corrections.  In  principle,  the  initial-guess  field  can  be  generated  in  a  variety  of  ways. 
An  accurate  forecast  of  the  field  being  analyzed  makes  the  most  ideal  initial  guess. 

Two  rather  reasonable  approaches  are  to  apply  a  surface-fitting  tcciinique  to  the  data 
or  to  use  a  recent  analysis  (e.g.,  two  hours  earlier),  if  one  is  available,  or  some 
combination  or  "blend”  of  the  two  approaches.  Our  version  of  the  program  uses  these 
surface  fitting  and  persistence  approaches. 

In  the  surface  fitting  procedure,  the  three-dimensional  surface  is  sought  which 
best  fits  the  observed  field  of  data.  The  equation  is  of  the  form 

<t>  =  a0  4  Vi  4  a2X2  +  a3X3  +  •  '  •  4  apxp  (3-3) 

where  the  dependent  variable  <p  is  the  variable  to  be  analyzed  over  an  evenly-spaced 
field  of  grid  points,  the  a’s  are  the  coefficients,  and  the  x’s  are  the  independent  vari¬ 
ables  which  are  functions  of  the  locations  of  the  observations  of  <p,  i.e.,  x,  y,  and  z. 
Because  there  are  a  large  number  of  terms  that  can  be  generated  from  the  possible 
combinations  of  x,  y,  and  z,  many  of  which  may  not  be  significant  in  explaining  the 
variability  of  4> ,  it  seems  reasonable  to  attempt  to  reduce  the  number  of  possible 
“predictors”  through  the  use  of  a  systematic,  stepwise  screening  procedure.  The 
method  chosen  is  that  of  screening  regression.  From  an  array  of  possible  predictors 
(x, y ,  z,  > y ,  xz,  etc.),  the  screening  procedure  first  selects  the  one  that  has  the  highest 
linear  ccrrelation  with  the  predictand  in  question.  This  predictor  is  then  held  con¬ 
stant,  and  partial-correlation  coefficients  between  the  predictand  and  each  of  the 
remaining  predictors  are  examined;  the  predictor  now  associated  with  the  highest 
coefficient  is  the  second  one  selected.  Additional  predictors  are  chosen  similarly 
until  a  selected  predictor  fails  to  explain  a  significant  additional  percentage  of  the 
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remaining  variance  of  the  predictand.  The  coefficients  a^,  a^,  etc.,  are  then  derived 
by  the  method  of  least  squares.  The  initial-guess  field  is  then  generated  by  solving 
the  derived  equation  at  each  of  the  grid  points  in  the  analysis  volume. 

The  program  is  written  so  that  it  is  possible  to  compute  an  initial  guess  from  a 
blend  of  persistence  and  surface  fitting.  This  is  done  to  incorporate  time  continuity, 
as  the  observations  are  taken  at  two-hour  intervals  on  alternate  days  between  0600  and 
1400  MST.  For  the  first  observation  time  of  the  day  (0600  MST),  however,  no  weight 
is  given  to  persistence  because  the  previous  available  observation  is  usually  forty  hours 
old.  The  formula  us  ed  in  the  program  to  calculate  the  initial-guess  value  of  <p  at  a  grid 
point  is  ' 

1  W 

0j  G  =  TTw —  ^(surface  fit)  +  ^(persistence)  (3-2) 

p  p 

where  W  =  weight  of  persistence.  By  assigning  a  value  of  zero  to  W  ,  persistence  is 
given  no  weight  and  the  initial  guess  is  based  solely  on  the  surface  fit.  A  value  of  one 
for  Wp  assigns  half  the  weight  to  each  while  a  very  large  value  of  essentially  gives 
all  the  weight  to  persistence. 

The  initial-guess  portion  of  the  computer  program  is  a  subroutine  called  IGUESS. 

3.3.3  Correction  Procedures  (subroutines  PCORR  and  TCORR) 

Because  the  arialysis  during  solution  may  only  be  defined  at  grid  points,  the 

information  from  observation  points,  which  in  general  are  at  random  locations  within 
grid  blocks,  must  be  translated  to  grid  points  to  permit  the  definition  of  internal  bound¬ 
ary  points.  For  each  observation,  the  difference  between  the  observed  value  and  the 
value  computed  for  that  location  by  interpolating  among  the  initial-guess  values  at  the 
surrounding  grid  points  is  computed.  This  difference  is  then  translated  to  the  nearest 
grid  point  and  added  to  that  grid  point  as  a  correction  to  the  initial  guess.  When  a  grid 
point  is  subject  to  multiple  corrections  (from  several  different  observations),  an  over¬ 
all  correction  is  computed  as  the  arithmetic  average  of  the  several  individual  correc¬ 
tions.  Corrected  grid  points  are  designated  as  internal  boundary  points  and  are  iden¬ 
tified  as  such  in  the  computer. 

3.3.4  Relaxation  Methods 

After  the  internal  boundary  points  have  been  defined,  the  next  step  in  CRAM  is  to 
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compute  interpolated  values  from  these  points  at  all  uncorrected  grid  points.  The 
interpolation  is  accomplished  by  requiring  that  the  values  at  the  uncorrected  grid 
points  satisfy  Poisson’s  equation: 

V2  0(i,j,k)  =  F(i,j,k)  (3-3) 

where  <p  is  the  value  of  the  analysis  parameter  at  grid  point  (i,j,k),  F  is  a  forcing 
function  defining  the  shape  of  the  field  of  0,  and  V  is  the  finite-difference  Laplacian 
operator  in  three  dimensions.  The  Laplacian  of  a  good  initial-guess  field  provides  a 
good  forcing  function.  In  the  program  (subroutine  FORCE),  the  formula  is  written  as 
follows: 

V2  0(i ,j ,k)  =  0(i  +  l,j,k)  +  0(i-l,j,k)  +  0(i,J  +  l,k)  +  0(i,j-l,k) 

-  40(i,j,k)  +  KA[0(i,J,k  +  l)  -  0(i,j,k)]  (3-4) 

-  Kg[0(i,j,k)  -  0(1, j,k-  1)) 


where  and  Kg  are  coefficients  pertaining  to  the  vertical  gradient  computed  above 
and  below  the  analysis  grid  point,  respectively. 

The  horizontal  and  vertical  space  gradients  have  been  expressed  separately  in 
this  equation.  There  are  several  practical  reasons  for  doing  this,  all  having  to  do  with 
the  scaling  of  the  vertical  coordinate  used  in  the  analysis  technique.  The  basic  hori¬ 
zontal  grid  length  is  10  km,  while  in  the  vertical  the  grid  length  is  about  an  order  of 
magnitude  smaller.  Measuring  gradients  over  this  smaller  distance  is  generally 
compensated  for  by  the  fact  that,  for  most  atmospheric  variables,  the  vertical  gradients 
are  much  larger  than  the  horizontal  gradients  so  that  one  may  be  justified  in  using 
Ka  =  Kg  =  1  in  Eq.  (3-4).  The  relative  inequalities  between  the  grid  distances  above 
and  below  the  grid  point  which  arise  from  the  uneven  spacing  of  the  artillery  zones 
may  be  accounted  for  by  assigning  different  values  to  K^  and  Kg.  Also,  it  is  possible 
in  the  program  to  consider  the  objective  analysis  problem  in  a  two-dimensional  frame¬ 
work  by  simply  setting  K^  =  Kg  =  0. 

The-  solution  of  Eq.  (3-3)  is  accomplished  by  relaxation.  Let  the  values  at  all  the 
boundary  points  (internal  and  perimeter)  remain  unchanged  throughout  the  relaxation 
computations.  For  each  non-boundary  point  in  the  grid  array,  the  residual,  R,  is 
computed  by 
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R(i- j»k)  =  V%(i,j,k)  -  F(i,j,k).  (3-5) 

New  estimates  of  <p  are  then  computed  by 

4>(  i.j.k)  =  aR(i,j,k)  +  $(i,j,k)  (3-6) 

where  a,  a  relaxation  coefficient,  is  an  input  parameter.  New  residuals  are  then 
computed.  The  iterative  procedure  continues  until  all  residuals  are  less  than  a 
specified  value  (e).  This  procedure  is  accomplished  in  subroutine  RELAX. 

3.3.5  Smoothing  (Subroutine  SMOOTH) 

A  smoothing  operator  has  been  developed  which  is  designed  to  eliminate  from 
the  analysis  undesirable  small-scale  features.  The  expression  used  is 
S[tf>(i.j.k)l  =  W^fi.J.k)  +  W2(*(i  +  l,j,k)  +  0(i  -  lj.k) 

+  *(i,j  +  l,k)  +  *(i,j-l,k)l  (3-7) 

+  W3{<Mi,j,k  +  l)  +  0(ifj,k- 1)1, 

/  ,  -  Bh  B 

w  = - 5 -  •  w  = - - - w  = - - - 

1  6  +  4B  +  2B  ’  2  6  +  4B.  +  2B  ’  3  6  +  4B.  +  2B 

hv  hv  h  v  / 

The  parameters  and  By  control  the  degree  of  smoothing  to  be  performed.  A 
value  of  zero  assigned  to  each  results  in  the  entire  weight  being  given  to  0(i,j,k),  i.e., 
no  smoothing.  Large  values  of  B^  and  By  result  in  heavy  smoothing.  Also,  the  smooth¬ 
ing  can  be  restricted  to  two  dimensions  by  setting  By  =  0. 

3-3.6  Ballistic  Computation 

The  CRAM  program  produces  analyses  of  u-  and  v-wind  components,  temperature, 
and  density  for  each  of  ten  artillery  zones.  These  zones  are  listed  in  Table  2-1  for 
an  analysis  area  17  grid  units  in  the  I-iirection  (from  1=  50  to  1  =  66)  and  15  grid  units 
in  the  J-direction  (from  J  =  345  to  J  =  359)  or  about  170  km  X 150  km.  A  portion  of  this 
grid  is  shown  in  Fig.  2-1.  The  analysis  area  is  purposely  made  large  so  that  analysis 
problems  associated  with  the  boundaries  are  far  enough  removed  to  not  contaminate 
the  analysis  in  the  region  of  interest,  i.e.,  near  the  firing  area.  The  ballistic  fields 
are  produced  within  the  program  by  applying  the  zone  weighting  factors  that  are  given 
in  Table  3-1.  A  minor  modification  in  the  program  will  permit  ballistic  computations 
of  "lines,,  other  than  merely  Line  10. 
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An  additional  feature  of  the  program  is  that  the  final  printed  output  for  the 
analysis  field  can  be  analyzed  with  contours  printed  for  any  desired  isopleth  interval. 

If  desired,  portions  of  the  analysis  cycle  may  be  repeated  several  times.  The 
generated  analyses  may  be  used  as  an  improved  initial-guess  field,  with  subsequent 
forcing,  correcting,  relaxing,  and  smoothing  to  make  the  analysis  fit  the  observations 
more  closely.  This  cycling  feature  can  be  employed  for  error  checking  purposes  by 
discarding  observations  which  differ  from  the  analysis  by  a  specified  value.  This  value 
may  be  large  for  the  first  cycle  with  successively  smaller  values  used  for  subsequent 
cycles . 

3.4  Production  of  Analyses  by  Application  of  CRAM 

The  various  options  and  parameters  discussed  above  were  given  a  limited 
amount  of  testing  on  a  portion  of  the  data  sample.  A  much  more  extensive  testing 
would  be  required  in  order  to  achieve  the  optimum  results.  Nevertheless,  within  the 
constraints  of  the  present  study,  it  was  possible  to  produce  meaningful  analyses.  In 
all,  there  were  80  observation  times  to  be  analyzed  (16  days,  5  observation  times  per 
day)  for  10  zones  for  four  variables  (u-  and  v-wind  components,  temperature,  and 
density),  making  a  grand  total  of  3,200  maps  that  were  analyzed  using  the  computer. 

The  initial-guess  procedure  of  blending  a  surface  fit  with  persistence  was 
employed,  using  W  =  1/2  in  Eq.(3-2).  The  forcing  function  (F)  was  set  equal  to  the 

p 

Laplacian  of  the  initial  guess. 

The  values  of  K  and  K  were  chosen  to  be  nearly  equal  to  one  (actually  varying 
A  B 

between  0.81  to  1.3  to  account  for  the  uneven  zone  spacing),  causing  the  analysis  to  be 

made  in  three  dimensions  rather  than  two.  The  number  of  iterations  required  was 

near  a  minimum  using  a  relaxation  coefficient  of  a  =  0.25.  (For  two  dimensions, 

a  -  0.22  should  be  used.)  The  smoothing  operator  was  employed  using  B,  =2.0  and 

h 

B^  =  2.0,  which  gives  one-third  of  the  weight  to  the  grid  point  being  smoothed  and 
two-thirds  to  the  surrounding  six  grid  points. 

Analyses  were  produced  for  each  of  the  ten  artillery  zones  using  the  procedures 
outlined  above.  These  zone  analyses  were  then  combined  into  ballistic  analyses  by 
computing  the  weighted  sums  using  the  weights  given  in  Table  3-1.  It  should  be 
emphasized  that  these  analyses  are  the  result  of  initial  investigations  and  that  further 
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work  in  this  area,  such  as  performing  the  force-relax-smooth  cycle  more  than  once, 
would  likely  lead  to  improved  analyses.  Figures  3-2  through  3-6  show  the  five  analyzed 
fields  of  ballistic  u-  and  v-wind  components  (in  knots)  for  0600,  0800,  1000,  1200,  and 
1400  MST.  The  temperature  and  density  analyses  (departure  from  Standard,  in  percent) 
are  given  in  Figures  3-7  through  3-11. 

The  u-component  ballistic  wind  analyses  show  a  gradual  increasing  trend  over 
the  southern  half  of  the  analysis  area,  with  only  slight  changes  with  time  in  the  northern 
portion.  This  causes  the  contours  that  are  originally  oriented  north-south  to  become 
nearly  west-east  by  1400  MST.  There  is  also  a  gradual  increase  in  the  v-component 
value. 

The  temperature  and  density  analyses  (Figs.  3-7  through  3-11)  show  relatively 
small  gradients,  with  the  range  of  values  over  an  entire  map  being  generally  less  than 
1%  of  Standard.  The  temperature  values  show  a  cooling  trend  through  the  series, 
reflecting  a  cooling  trend  in  the  upper  zones.  The  surface  heating  affecting  the  lower 
zones  does  not  contribute  much  to  the  overall  ballistic  temperature  because  of  the  way 
the  individual  zones  are  weighted  (see  Table  3-1).  There  is  practically  no  spatial 
variability  in  density  until  the  final  map. 
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Fig.  3-3.  Ballistic  wind  analysis  using  CRAM:  0800  MST,  20  January  1965 


ME 


SIS  US 


ary 


,  20  January  1965  (percent 


Fig.  3-11.  Ballistic  temperature  and  density  analysis  using  CRAM;  1400  MST,  20  January  1965  (percent 


4.0  CALCULATION  OF  RESIDUAL  ERRORS 


The  main  purpose  of  developing  an  objective  analysis  technique  for  this  study  is 
to  produce  a  representative  estimate  of  the  meteorological  environment  through  which 
an  artillery  projectile  travels.  It  is  hoped  that  the  analysis  technique,  which  incor¬ 
porates  the  map-wind  concept,  will  be  belter  than  the  ballistic  winds  derived  from  a 
single  sounding.  One  meaningful  way  to  evaluate  the  relative  merits  of  two  or  more 
methods  of  computing  ballistic  winds  is  to  make  a  comparison  of  residual  errors  using 
actual  artillery  firings. 

Applying  corrections  to  the  firing  data  supplied  for  this  study  required  the  use 
of  unit  effect  data  supplied  by  BRL  (Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground)  based  on  firing  table  information.  Ballisticians  have  derived  these  unit 
effects  theoretically  by  first  computing  a  trajectory  from  equations  of  motion  using 
standard  conditions.  They  then  vary  each  of  the  variables  from  standard, one  at  a  time, 
by  a  unit  amount  (1%  for  density  or  temperature,  one  knot  for  range  wind  or  cross 
wind)  and  note  the  resulting  difference  in  range  or  deflection.  These  differences  are 
the  unit  effects. 

For  this  study  it  was  assumed  that  corrections  for  muzzle  velocity,  projectile 
weight,  powder  temperature,  etc.,  were  already  accounted  for  in  the  computation  of 
the  range.  Because  the  firings  were  conducted  at  fixed  azimuths  and  quadrant  ele¬ 
vations,  the  only  corrections  necessary  to  "move”  the  observed  impact  points  in  to 
the  target  were  temperature,  density,  and  range  w'ind  for  the  range  correction,  and 
drift,  coriolis  force,  and  cross  wind  for  the  deflection  correction.  The  drift  and 
coriolis  force  are  independent  of  the  meteorological  values  used,  so  that  only  ballistic 
wind,  temperature,  and  density  are  pertinent  in  drawing  comparisons  among  various 
analysis  techniques. 

4.1  Residual  Errors  Using  CRAM  and  Ft.  Huachuca  Sounding 

The  provided  gun  data  and  unit  effects  were  used  in  computing  artillery  correc¬ 
tions.  (Only  results  using  Gun  No.  1  are  presented  here  because  the  characteristics 
of  Gun  No.  2  were  very  similar.)  Using  CRAM  analyses  concurrent  with  the  24  firing 
times,  values  of  ballistic  temperature,  density,  and  wind  (resolved  into  range  and 
cross  wind  components)  valid  at  the  midpoint  of  the  firing  trajectory  were  extracted. 
The  results  are  shown  in  Table  4-1,  where  the  odd-numbered  cases  refer  to  Gun  No.  1. 
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TABLE  4-1 

ARTILLERY  CORRECTIONS  BASED  ON 
CRAM  ANALYSIS 
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TABLE  4-1  (continued) 
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The  target  (hypothetical)  coordinates  are  given  in  metors  in  the  Universal  Transverse 
Mercator  grid  system.  The  range,  deflection,  and  vector  errors  designated  “NO 
CORRECTIONS”  arc  based  on  the  distance  between  the  target  and  the  given  center  of 
impact  (Cl).  The  correction  due  to  any  one  meteorological  parameter  is  obtained  by 
multiplying  its  value  by  the  unit  effect  for  that  parameter.  The  range,  deflection,  and 
vector  errors  "WITH  CORRECTIONS”  are  residual  errors  based  on  the  corrected 
impact  point  after  the  drift,  coriolis,  and  meteorological  effects  have  been  accounted 
for.  (Negative  range  errors  refer  to  impact  short  of  the  target  and  negative  deflection 
errors  refer  to  impact  to  the  leit  of  the  target  when  viewed  from  the  gun  location.) 

Also  given  are  the  x  and  y  components  of  the  vector  error  which,  when  added  to  the 
target  coordinates,  give  the  coordinates  of  the  corrected  impact  point.  The  corrections 
due  to  temperature,  density,  and  range  wind  affect  only  the  range  (negative  corrections 
decrease  the  range),  while  corrections  due  to  cuss  wind  affect  only  the  deflection 
(negative  corrections  displace  the  impact  point  to  the  left).  The  variations  in  target 
coordinates  and  unit  elfects  are  due  to  adjustments  made  for  non-meteorological 
factors  by  BRL. 

Ballistic  quantities  determined  using  the  Ft.  Huachuca  sounding  alone  were  also 
applied  to  the  firing  data.  Of  the  24  firing  times,  there  were  5  in  which  the  data  were 
not  available  to  construct  a  Line-10  ballistic  message.  The  results  of  the  19  cases 
are  shown  in  Table  4-2. 

4.2  Compa rison  of  Re sidual  Errors 

A  tabulation  of  average  range  and  deflection  errors  for  the  CRAM  and  Ft.  Huachuca 
data  showed  rather  conclusively  a  systematic  error  in  range  of  about  190  m,  the  mean 
Cl  being  about  190  m beyond  the  target.  Discussions  with  ballistics  experts  at  BRL 
were  held  to  see  if  the  reason  for  the  error  could  be  detected  and  corrected.  It  was 
brought  out  that  the  ballistic  coefficient  used  in  the  solution  of  the  equations  of  motion 
was  not  accurately  known  for  high  angle  firing  situations,  i.e.,  quadrant  elevation  above 
45  degrees.  One  reason  for  this  is  that  when  the  projectile  is  oriented  at  a  high  angle 
as  it  leaves  the  gun  barrel,  it  is  subjected  to  a  “summital  yaw”  in  the  vicinity  of  the 
top  of  the  trajectory.  This  results  in  an  erratic  descent  to  the  ground  and,  thus,  an 
unreliable  trajectory.  The  ballistic  coefficient  when  determined  by  a  relatively  small 
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number  of  high  nrgb'  fi  pint's  is  therefore  subject  to  error.  Because  the  quadrant 
elevation  for  the  Ft.  Hunt  i  c'ca  firings  was  of  the  order  of  G.T,  a  suggestion  was  made 
to  have  the  ballistic  coefficient  computed  b>  BRL  using  the  Ft.  Huachuca  firing  and 
meteorological  data.  The  detailed  results  were  not  available  in  time  for  publication, 
but  it  was  learned  that,  based  on  about  one  half  of  the  data  sample,  there  were  discrep¬ 
ancies  in  the  ballistic  coefficient  of  about  5  to  10%.  This  was  not  enough  to  change  the 
unit  effects  significantly,  but  it  did  cause  changes  in  range  (target  location)  of  as  much 
as  000  m. 

Rather  than  complete  elimination  of  any  comparison  of  range  residual  errors,  it 
was  decided  ns  an  interim  measure  to  remove  some  of  the  systematic  error  empirically 
A  linear  regression  equation  relating  range  error  to  range  wind  was  derived  by  the 
method  of  least  squares  on  the  combined  sample  of  24  CRAM  cases  and  the  19  Ft. 
Huachuca  cases.  The  data  are  plotted  in  Fig.  4-1.  The  correlation  coefficient  was  -0.82 
The  equation  for  the  regression  line  is 

C  .122.87  -  7.11  (RW)  (4-1) 

where  C  is  the  correction  factor,  in  meters,  to  be  subtracted  from  the  conventionally 
determiner!  range  error  and  RW  is  the  range  wind  in  knots.  Corrections  made  by  this 
method  are  in  good  qualitative  agreement  with  those  made  to  date  by  recomputing  the 
ballistic  coefficient. 

Table  4-1  shows  a  summary  of  residual  errors  for  CRAM  and  the  single  sounding 
from  Ft.  Huachuca  where  the  range  residual  error  has  been  empirically  corrected.  The 
average  absolute  errors  are  based  on  the  19  cases  where  comparison  is  possible.  The 
average  range  and  deflection  errors  are  9.G  m  and  9.9  m  smaller  for  CRAM  than  for 
Ft.  Huachuca.  The  residual  vector  difference  is  12.1  m  lower  for  CRAM.  The  average 
residual  range  errors  before  the  empir  cal  corrections  were  231.1  m  for  CRAM  and 
243.6  m  for  Ft.  Huachuca.  The  differences  in  deflection  errors  are  significant  at  the 
15%  level  when  a  Student's  T-test  for  paired  comparisons  (Wadsworth  and  Bryan,  1963) 
is  employed. 

The  residual  errors  are  shown  graphically  in  Fig.  4-2. 

Differences  in  residual  errors  between  the  two  methods  should  be  the  largest 
when  winds  are  strong  and/or  considerable  horizontal  gradient  in  the  meteorological 
parameters  exists. 
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ARTILLERY  CORRECTIONS  BASED  ON  FT.  H 
METEOROLOGICAL  SOUNDINGS* 
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TABLE  4-2  (continued) 
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TABLE  4-2  (continued) 
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TABLE  4-2  (continued) 
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TABLE  4-2  (continued) 


:'A 


•  o 

a  • 

1  >  5  R 

X 

•  o  n 

.  x 

(  o  — 

X 

i 

X 

i 

X 

i  •  r 

X 

i  a  • 

X 

in  x  x  in 

X 

z  o  o 

X 

O  u  CM 

X 

•  ■* 

X 

R 

X 

u  z  a  ax 

X 

a  r  o  • 

z 

X 

a  a  a  r 

o 

X 

a  a  a  r 

—  ?  a  <r  3 

3 

X 

o  a  a  i 

r  a  3  n 

• 

X 

u  a 

u  *  •  •  • 

G 

X 

Z 

a  n  o  o  aj 

X 

o 

a  a  3 

a  i  g  n  a 

R 

X 

X  o  o  • 

a  r  aj 

X 

H 

R  2  3  <r 

O  1  1 

X 

< 

—  <  a  n 

u 

II 

X 

h- 

3  a  a  i\i 

Z 

X 

i 

o 

X 

o 

1 

i  •  cr  aj 

R 

X 

X 

UJ 

1  I-  o  • 

u 

X 

X 

i  u  a  r 

R 

a 

X 

***** 

i  a  a  <r 

U 

a 

X 

>  a  aj 

a 

a 

X 

a  n  o  n  g 

o 

X 

1/3 

a  a  —  aj  on 

u 

X 

UJ 

z  a  n 

a  •  •  •  • 

X 

z 

R  O  • 

—  -  c  « 

R 

X 

o 

•  a  a  n 

R  — •  G  Aj  — • 

L 

X 

N 

<n  a  ~  o 

-  1  II 

X 

z  a  a  a 

Z 

a 

X 

o 

O  Q 

3 

o 

X 

— 1 

X 

R  a  a  s 

X 

X 

o  e  o  • 

X 

o 

a  z  a  x 

X 

<t 

a  •!  3  o 

a 

X 

D 

g  ~  ■  ■  ■*' 

3 

X 

X 

o  “  “ 

a 

X 

u 

X 

•3D 

>  n  a 

X 

o>-o* 

nan  — 

X 

z  u  a  3 

•  •  f  •  . 

3 

X 

•  a  a  - 

a  go 

<J 

X 

>  a  - 

O  —  AJ 

. 

X 

— 

a 

C 

X 

R 

n 

X 

a 

X 

V 

X 

a 

X 

— 

ii 

X 

• 

z 

X 

AJ 

o 

X 

UJ 

a 

X 

{/) 

on  x 

a 

R 

X 

< 

n 

a  a 

a 

u 

X 

u 

C\i 

R  G 

a  a 

a 

X 

R  <  3 

S3  GO 

3 

X 

2 

a  z 

<3  Z  Z 

a 

X 

J 

o  - 

a  r  — 

o 

X 

o 

a  q  z 

<  <  >  ?  3 

u 

X 

a 

<  a  o 

a  a  r 

X 

R  o  • 

a  -  a  on 

on 

X 

lO 

o  o 

•  a  or  o  3 

X 

tii 

vO 

U  R 

a  5  z  z  o 

j 

X 

K 

2 

3 

o  a  a  <  a 

o 

X 

< 

< 

R 

a  i-  n  a  u 

X 

Q 

~> 

a 

R 

a 

X 

h* 

a 

o 

X 

rj 

n 

5 

u 

X 

•  in 

a  • 

>50 
o  a 
u  - 


i  a  • 

on  X  S  3 

z  o  nj 

o  u  m 


o  z  a 

R 

a  r  o 

• 

z 

a  j  a 

AJ 

3 

a  a  a 

AJ 

— 

o  a  a 

— 

R 

u  Q 

1 

o 

z 

a 

o 

a  a 

a 

a 

— 

IOC 

• 

a 

R 

r  z  a 

3 

o 

< 

-  <  a 

G 

u 

R 

3  a  a 

AJ 

on 

i 

G 

| 

i  •  a 

G 

a 

1  R  o 

• 

X 

i  o  a 

3 

R 

•— 

i  a  a 

a 

U 

a 

>  a 

AJ 

a 

on 

a 

a 

a 

z  a 

X 

a 

z 

R  O 

• 

o 

i  j  a 

<T 

R 

N 

m  a  a 

<T 

— 

z  a  a 

r* 

“X 

c 

o  o 

3 

— 

•— > 

r  a  a 

O 

I 

u  o  o 

• 

J 

uj  z  a 

— 

< 

a  <t  a 

a 

3 

a  a  a 

3 

3 

r 

o 

a 

0 

< 

•  a 

R 

> 

o  R  o 

• 

z  u  a 

3 

• 

i  a  a 

3 

a 

>  a 

C 

o 

_ 

a 

R 

a 

~  g  G 

R  i>  3 


X)  n 
a  a 


in  m  o 
n  nj 

i  t 


-  on  - 
o  -  n 

•  •  •  • 

<M  O  O  <T 

-  g  m  - 

i  i  i 


n  ?  n  - 

•  •  •  * 

>0  G 

-  AJ 


U 

o 


UJ 

<  <j 

U  aj 

Z 

3  nj 
O 

a 


n 

o 

z 

< 

o 

r 

M 


CD 

a 

n 

C 

G 

XI 


o  z 
cr  g 

o  • 

O  G 
U  CD 
3 
R 
O' 

n 


a 

UJ 

R 

Lu  uJ 

i  a 
< 
a 
< 
a 


a 

o 

UJ 

t- 

UJ 

5 


a  a 
z  z 

>33 


-  UJ  G 
l/l  O  Ul 

z  z 
a  a  <  a 
r  a  a  u 


o 


xi 

XI 

<r 

r 


50 


CORIOLIS  CORRECTION-  30.02  ORIFT  CORRECTlON= 


TABLE  4-2  (continued) 
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TABLE  4-2  (continued) 
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CRAM  analysis 


Ft.  Huachuca 


ate/time 

Case 

Vector 

error 

Range* 

error 

Deflection 

error 

Vector 

error 

Range* 

error 

Deflection 

error 

/ 1 6 

AM 

1 

28.9 

-  18.8 

22.0 

45.5 

-  3.1 

PM 
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-  68.7 

-  60.8 

122.4 

-  68.9 

/18 

AM 

5 

47.0 

12.5 

-  45.3 

75.8 

11.9 
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PM 
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20.1 

12.0 

-  16.1 
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AM 

9 

112.6 

112.6 
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45.2 

PM 
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390.0 
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387.0 
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503.2 

381.5 
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— 

PM 
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‘With  empirical  correction. 
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Deflection  error 


Fig.  4-2.  Residual  errors  for  artillery  firings  (Gun  No.  1)  based  on  CRAM 
and  Ft.  Huachuca  (from  Table  4-3). 
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5.0  TIME  VARIATIONS  IN  THE  BALLISTIC  WIND 


5.1  Computation  Procedure 

Variability  measures  appropriate  to  the  ballistic  wind  data  and  to  the  way  in 
which  the  basic  rawinsonde  observations  were  scheduled  are  the  root-mean-square 
(rms)  wind  component  differences  and  the  rms  difference  vector  magnitudes. 

These  quantities  were  computed  at  all  grid  points  for  the  time  lags  of  2,  4,  6, 
and  8  hours,  and  for  lines  2,  4,  6,  8,  and  10.  The  sample  size  varied  with  time  lag 
because  of  the  procedure  used  in  scheduling  the  observation  runs.  Table  5-1  shows 
the  distribution  of  differences  computed  by  lag  and  terminal  hour,  for  any  given  line. 

TABLE  5-1 

DISTRIBUTION  OF  WIND  DIFFERENCES  COMPUTED  BY  LAG  AND  TERMINAL  HOUR 


Lag 

(hr) 

Terminal  hour  (MST) 

Row 

totals 

0800 

1000 

1200 

1400 

2 

16 

16 

16 

16 

64 

4 

X 

16 

16 

16 

48 

6 

X 

X 

16 

16 

32 

8 

X 

X 

X 

16 

16 

The  computer  output  consisted  of  a  field  of  grid-point  values  of  the  variability 
measures  for  each  terminal  hour  and  lag,  and  pooled  values  for  all  terminal  hours  for 
each  lag.  Subsequently,  field-averages  were  calculated  from  the  computer  output. 

These  field-averages  were  obtained  from  that  portion  of  the  analysis  field  in  which 
the  data  are  densest,  extending  from  grid-point  347  through  355  inclusive  in  the  north- 
south  direction,  and  from  grid-point  52  through  64  inclusive  in  the  east-west  direction. 
Details  of  the  computational  procedure  are  given  in  the  computer  program  specifications 
that  follow. 

5.1.1  Computations  for  Time  Variability  Studies 

These  specifications  outline  the  computation  of  the  rms  vector  difference  between 
pairs  of  ballistic  winds  .-.eparatedby  given  time  intervals,  for  a  set  of  time  intervals. 

The  statistics  are  computed  for  selected  line  ballistic  winds  and  at  each  gnd  point  in 
the  analysis  fie.d.  Also  computed  are  various  quantities  to  be  used  in  predictability 
studies. 
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5.1.2  Input 

The  following  quantity  is  input: 

l  J  (L,  i,  j,  t)  The  two-component  vector  ballistic  winds  (knots)  for  the  line  L. 

at  the  time  t,  and  the  horizontal  grid  point  with  coordinate  i.j, 
where 

L  =  2,  4,  6,  8,  10 

t  =  1/12-0600  MST .  2/12-1400  MST 

i  =  1 . I 

j  =  1 . J- 


5.1.3  Output 

The  following  quantitie 

a  (t.  L,  i.  j,  q)  Root-mean-square  ballistic  wind  component  difference  (knots) 
over  the  time  interval  r,  ending  at  the  observation  hour  q, 
where, 

r  -  2,  4,  6,  8  hours 
q  =  1  +  t/2  (r  r  2,4 .6.8) 

L,  i,  j,  q)  Same  as  immediately  preceding,  but  for  v-component. 

ct(7,  L,  i,  j,  q)  Same  as  immediately  preceding,  but  for  vector  difference. 

a  (r.  L.  i,  j)  Same  as  cr  above,  but  for  all  observation  hours, 
u  u 

av(~>  L,  i,  j)  Same  as  <7  above,  but  for  all  observation  hours. 
ct(t,  L,  i,  j)  Same  las  a  above,  but  for  all  observation  hours. 

The  following  quantities  are  to  be  output  with  the  predictability  experiment 
output  (see  Section  8.0): 


o  (r,  L,  i,  j,  q)  Root-mean-square  u-component  error  of  selected  persistence 
U,f>  forecasts  over  time  period  r,  and  valid  at  the  observation 

hour  q. 

q  =  2  +  r/2  (t  =  2.4,6) 


(t,  L,  i,  j,  q)  Same  as  immediately  preceding,  but  for  v-component  error, 

u  (r,  L,  i,  j,  q)  Same  as  immediately  preceding,  but  for  vector  error. 

a  (t,  L,  i,  j)  Same  as  o  above,  but  for  all  valid  times. 

u, p'  u,p 

a  (t,  L,  i,  j)  Same  as  cr  above,  but  for  all  valid  times. 

v , p  v,p 

°p(r.  L,  i.  j)  Same  as  cr^  above,  but  for  all  valid  times. 
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The  following  quantities  are  to  be  used  in  the  predictability  experiments: 

D^fT,  I,,  i,  j,  q)  u-component  differences 
D^(r,  L,  i,  j,  q)  v-component  differences 
q  -  1  +  t/2  (t  =  2, 4. 6, 8) 

5.1.4  Computations 

The  sequence  of  map  times  is  to  be  indexed  (p,q)  where  the  index  p  represents 
the  observation  day,  p  =  1,  ....  16,  and  the  index  q  represents  the  observation 
hour,  q  =  1,  2,  5,  4,  5. 

Then,  for  fixed  L,  p,  i,  j,  the  difference  components 
Du(L,  p,  i.J,  T,q)  •  Uq  -  Uq_r/2 

Dv(L.  p,  i,j,  r,  q)  =  Vq  -  Vq_r/2 

[(1  <  7/2)  £  q  ^  5] 
are  computed. 

For  each  q,  the  mean  square  differences  over  all  p  are  computed: 

2  1^2 
uu(7,  L.  i,hq)  =  L  Du(r,  L,  i,  j,  p,  q) 

P=1 

1  f6  2 

CTV.(7.  L,  i,  j,  q)  =  “  ^  Dv(r,  L,  i,  j,  p,  q) 

P=  1 

2  2  2 
a  (r,  L,  i,  j.  q)  -  ^ 

[(1  >  p  -  q  <  5] 

and  the  mean  square  errors  of  the  selected  persistence  forecasts, 

2  ,  ...  . 

CTu  n  ~ ’  L’  J’  J’  q) 

2  .  ...  . 
a  (7 ,  L,  i,  j,  q) 
v,p 

a  (7,  L,  l,  j,  q)  , 

P 

are  obtained  by  the  same  formulas,  but  with  q  restricted  to 
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The  corresponding  root-mean-square  output  quantities  are  obtained  by 

1/2 

9 

cr  -  (a  ) 


Now,  for  fixed  ",  L,  i,  j,  the  mean  square  differences  over  all  hours  q  are  obtained 
from  the  above  by 


a“(r,  L,  i,  j)  - 


2  2 

and  similarly  for  a  and  cr  . 

u  v 


1  r  2, 

— r  Y  °  <7 

5-  .  t 


,  L,  i,  j,  q) 


q=1+T 


The  mean  square  persistence  forecast  errors  are  obtained  by 


L.  i.  j) 


2  2 

and  similarly  for  tr“  and  cr 

u.p  v,p 


2 


a 


(r,  L,  i,  j,  q) 


5_. 2 _ Ti me-variability  as  a  Function  of  Lag 

5.2.1  Variability  of  Line- 1 0  Ballistic  Winds 

Line-1 0  component  variability  parameters  obtained  from  objective  analyses  at 
the  midpoint  of  the  test-firing  trajectory  are  shown  in  Figs.  5-1  and  5-2.  For  these 
firings  the  u-component  corresponds  approximately  to  the  range-wind  component  and 
the  v-ccmponent  to  the  cross-wind  component.  As  may  be  seen  from  Table  5-1,  there 
is  a  systematic  reduction  of  sample  size  with  time  lag,  which  would  tend  to  produce 
progressively  more  uncertainty  in  the  variability  estimate  as  time  lag  increases. 

A  qualitative  estimate  of  the  effect  of  the  sampling  inconsistency  on  the  slope  of 
the  variability-lag  curve  was  obtained  by  reducing  the  sample  size  for  the  2-,  4-,  and 
6-hr  lag  variability  estimates  to  that  of  the  8-hr  lag  estimates.  This  was  done  by 
including  only  differences  between  the  1400  ballistic  wind  and  other  ballistic  winds  in 
the  sample.  The  reduced  sample  variability  is  also  shown  in  Figs.  5-1  and  5-2.  When 
compared  to  the  full  sample,  a  slight  but  systematic  decrease  in  variability  is  evident 
which  is  more  pronounced  in  the  u-component  variability. 
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— -  Full  sample 

•  reduced  sample 

•  Huachuca  single-station  wind 

J _ I _ J _ 

4  6  8 

Lag,  hr 


time  variability,  v -component . 


( 


The  variability  parameters  obtained  for  the  single  sounding  at  Ft.  Huachuca  and 
the  full  sample  are  also  shown.  As  would  be  expected,  the  map  ballistic  winds  show 
smaller  rms  differences  than  do  the  single-sounding  winds.  This  would  imply  a  slower 
rate  of  decay  of  map  wind  with  time,  for  artillery  correction  purposes.  However,  the 
variability  parameter  values  exceed  5  knots  in  both  components  by  the  4th  hour,  for 
the  lfi-day  sample.  Approximate  unit  effects  applicable  to  the  test  firings  would  give 
rms  displacement  errors  of  70  m  and  100  m  in  deflection  and  range,  respectively,  for 
5  knot  rms  wind  differences. 

5.2.2  Comparison  with  Ft.  Sill  Experiment  Results 

Figures  5-3  and  5-4  show  the  variability  as  a  function  of  lag  for  the  Line  4 
ballistic  wind  components  for  an  8-day  sample  at  Ft.  Sill,  as  obtained  from  Bellucci 
(19 61).  Also  show;  are  Line  4  rms  component  differences  for  the  data  obtained  in  the 
present  study.  In  these  figures  the  field-average  variability  (sec  Section  5.1)  is  shown 
for  the  analysis  region,  as  well  as  the  variability  in  selected  single  soundings.  The 
lower  variability  characteristic  of  space-averaged  map  ballistic  winds  as  compared  to 
single  sounding  ballistic  winds  is  evident  in  both  figures.  Terrain  effects  on  the  varia¬ 
bility  of  low  level  ballistic  winds,  which  will  be  discussed  in  more  detail  in  Section  5.3 
below,  are  seen  in  the  generally  higher  u-component  variability  (Fig.  5-3)  at  Ft.  Huachuca 
than  at  Ft.  Sill.  The  v-component  variability  at  Douglas  is  noticeably  lower  than  at 
Ft.  Sill  in  Fig.  5-4.  The  significance  of  these  differences,  of  course,  is  questionable 
because  of  the  limited  sample  size  obtained  in  the  two  experiments.  However,  the 
irregularity  in  the  FI.  Sill  curve  does  suggest  an  underlying  phenomenon  characteristic 
of  the  western  Great  Plains  which  would  be  absent  in  mountainous  terrain.  Thic  is  the 
breakdown  of  the  low-level  nocturnal  jet  during  the  morning  hours.  The  Ft.  Gill  van- 
ability  at  4-hrlag  is  computed  from  differences  between  the  0600-  and  1000-CST 
ballistic  winds,  as  well  as  from  differences  between  0800-  and  1200-CST  ballistic  winds. 
Where  the  nocturnal  jet  is  present,  these  periods  would  be  those  of  relatively  large 
change  in  the  morning  hours  (to  which  the  observational  sample  was  confined). 

The  curve  for  the  relationship 
1/2 

a  ot  '  , 
t 

where  o  is  in  knots  and  t  in  hours  (Bellucci,  1961), is  also  shown  in  Figs.  5-3  and  5-4. 
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Fig.  5-3.  Line-4  time  variability,  u-componcnt. 
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This  relationship  appears  more  consistent  with  the  mountainous  terrain  rms  wind 
component  differences  than  with  Bellucci’s  variability  for  the  Ft.  Sill  data,  and  agrees 
better  with  the  field-average  map-wind  variability  than  with  the  single-sounding  vari¬ 
ability.  Note,  however,  that  this  relationship  does  not  fit  the  Line-lOdata  well  (see 
Figs.  5-1  and  5-2). 

5.3  Spatial  Distribution  of  Time-variability  Parameters 

5.3.1  Effects  of  Analysis  Procedure  on  Computed  rms  Differences 

The  method  used  for  obtaining  the  initial-guess  zone  wind  fields  for  the  analysis 
procedure  obviously  has  influence  on  the  computed  fields  of  time-variability  parameters. 
Use  of  the  preceding  map  for  the  initial-guess  field  would  tend  to  produce  maximum 
rms  differences  at  the  data  points,  while  use  of  a  second-order  surface  of  best  fit 
would  tend  to  produce  maximum  variability  at  those  boundary  points  furthest  from 
dato-dense  regions. 

Initial-guess  fields  were  actually  generator  by  combining  the  2-hr  previous  map, 
where  one  was  available,  with  the  best  fit  surface,  with  equal  weights.  Thus,  as  would 
be  expected,  the  two  features  described  above  are  evident  in  the  grid-point  fields  of 
rms  differences. 

These  features  are  most  pronounced  for  the  lower  level  ballistic  winds,  where 
the  data  points  are  clustered  about  the  fixed  balloon  release  points.  Figures  5-5  and 
5-6  show  typical  fields  of  rms  vector  difference  magnitude,  in  this  case  for  time  lag 
2  hr,  as  obtained  for  Line  2  and  Line  10.  Note  the  contrast  between  the  two  fields 
caused  by  the  more  widespread  distribution  of  data  points  used  to  obtain  the  Line- 10 
analyses . 

The  local  maxima  evident  in  the  low-level  fields  tend  to  become  m»  re  pronounced 
with  increasing  time  lag.  Figure  5-7  shows  the  field  of  vector  variability  for  Line  2 
based  on  an  8-hr  time  lag,  which  may  be  compared  with  Fig.  5-5. 

5.3 .2 _ Horizontal  Distribution  of  Time  Variability  Parameters 

Figure  5-6  is  typical  of  the  upper  level  rms  ballistic  wind  difference  fields  by 
being  remarkably  uniform  over  most  of  the  analysis  area.  In  fact,  the  range  in  vari¬ 
ability  over  the  data-dense  portion  of  the  field  specified  in  Section  5.1.1  is  approxi¬ 
mately  one-half  knot  at  Line  10,  for  2-hr  time  lag,  and  increases  to  about  one  and  a  half 
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knots  for  8-hr  time  lag. 

At  the  lower  levels,  however,  some  terrain  effects  on  the  time  variability  may 
be  noted.  If  only  those  grid-point  values  nearest  the  balloon  release  points  are  plotted, 
and  isopleths  of  equal  rms  wind  (component  or  vector)  differences  are  drawn  with  a 
topographic  map  as  background,  maps  similar  to  those  illustrated  in  Figs.  5-8,  5-9, 
and  5-10  may  be  obtained  for  all  four  values  of  time  lag.  The  maps  shown  are  for 
8-hr  time  lag,  and  thus,  show  generally  larger  magnitude  rms  differences  than  maps 
for  shorter  time  lags  would  exhibit.  However,  the  patterns  evident  are  characteristic 
of  the  shorter  time  lags  also.  A  region  of  maximum  variability  on  the  v-component 
is  seen  in  the  triangle  Hereford,  Tombstone,  Fairbanks,  near  the  head  of  a  large 
north-south  valley  (see  Fig.  5-9),  and  lower  values  are  seen  in  the  Nogales-Patagonia- 
Sonoita  region  where  topographic  features  have  a  more  east-west  orientation.  The 
opposite  arrangement  of  maxima  and  minima  may  be  seen  in  the  u-component  vari¬ 
ability,  while  both  regions  exhibit  maxima  in  the  vector  wind  variability. 

The  range  in  variability  over  the  data-dense  portion  of  the  analysis  field  is 
generally  larger  at  Line  2  than  at  Line  10.  The  spatial  range  in  variability  is  about 
1  knot  at  Line  2  fcr  2-hr  time  lag.  and  increases  to  almost  2  knots  for  an  8-hr  time 
lag  (see  Fig.  5-10). 

5.3.3  Vertical  Distribution  of  Time  Variability  Parameters 

Vertical  profiles  of  the  time  var. ability  show  the  expected  general  increase  with 
height.  Figure  5-11  shows  such  profiles  for  field-average  rms  vector  difference 
magnitudes,  for  lags  of  2,  4,  6,  and  8  hr,  as  a  function  of  height.  Topographic  effects 
are  noticeable,  particularly  with  regard  to  short  period  variability,  in  the  irregularity 
of  the  rate  of  increase  with  height  in  the  lower  5000  m  in  Fig.  5- 1 1 .  Figure  5-12  shows 
similar  profiles  for  rms  u-component  differences,  in  which  the  topographic  effects 
on  the  low  level  profiles  are  less  pronounced.  Figure  5-14  shows  the  v-corw  onent 
differences  in  which  the  topographic  effects  are  most  evident. 

Some  evidence  for  the  association  of  these  low  level  effects  with  the  topography 
is  evident  in  Figs.  5-13  and  5-15,  in  which  selected  individual  grid-point  values, 
rather  than  field  averages,  are  shown.  Two  pairs  of  grid  points  closest  to  stations 
with  (a)  large  contrast  in  component  variability  as  indicated  by  the  Figs.  5-8  and  5-9, 
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and  (b)  a  similar  number  of  total  observations  in  each  sample,  were  selected  to 
illustrate  the  range  in  profiles  within  the  sample. 

Figure  5-13  shows  that  the  lower  portions  of  the  profiles  of  u-component  short- 
period  variability  represent  larger  variability  in  this  component  at  Huachuca  than  at 
Fairbanks.  Figure  5-15  shows  a  similar  contrast  in  the  v-component  profiles  for 
Hereford  and  Patagonia.  Note  the  decrease  with  height  in  v-component  variability  at 
Hereford  at  all  time  1---Ts  in  the  height  interval  between  Line  2  and  Line  4. 

5 .4  llour-to-hour  Range  of  Time  Variability  Parameters 

Figure  5-16  shows  the  range  of  variation  of  2-hr  and  4-hr  lag  variability  at  the 
grid  point  nearest  Ft.  Huachuca  at  Line  10  as  a  function  of  the  hour  of  the  day.  The 
hour  here  is  that  of  the  later  map  ballistic  wind  used  in  computing  the  vector  ballistic 
wind  differences. 

It  is  evident  that  the  range  in  variability  for  a  given  lag  is  definitely  smaller  than 
the  difference  in  variability  between  lags,  being  approximately  one-half  to  one  knot  as 
compared  to  four  knots.  Ar.  interesting  feature  to  be  investigated  in  later  studies  is 
the  peak  in  both  2-hr  and  4-hr  lag  variability  in  the  observations  made  at  1200  MST 
for  Line  10. 

This  feature  is  also  evident  in  the  Line- 2  variability-  for  2-hr  lag  at  the  same 
grid  point. 
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Field-  averaged  time  variability  (vector),  vertical  profile. 
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Fig.  5-12.  Field-averaged  time  variability  (u-componect),  vertical  profile . 
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6.0  SPATIAL  VARIABILITY  OF  BALLISTIC  WINDS 

To  examine  the  representativeness  of  ballistic  winds  estimated  from  observations 
somewhat  removed  from  the  environment  the  projectile  actually  passes  through, differ¬ 
ences  between  the  ballistic  wind  components  at  the  trajectory  midpoint  and  all  grid 
points  in  the  analysis  area  were  computed. 

Summary  maps  were  derived  for  all  80  sets  of  observations.  The  summary 
statistic  employed  is  the  rms  difference  of  the  ballistic  wind  components  at  the  tra¬ 
jectory  midpoint  and  all  grid  points  in  the  analysis  area.  For  example,  at  grid  point 
i,j,  the  rms  difference  for  the  u-component  of  the  ballistic  wind  is 


du«>-[mE<VV2. 


(6-1) 


where  UQ  is  the  u-component  of  the  ballistic  wind  at  the  trajectory  midpoint  and  U  is 
the  u-component  of  the  ballistic  wind  at  grid  point  i,j. 


Maps  of  the  rms  difference  for  the  two  components  (as  well  as  for  the  vector)  of 
the  ballistic  winds  were  prepared  for  Lines  2,  4,  6,  8,  and  10. 


These  maps  were  examined  along  with  observational  data  to: 


(a)  obtain  estimates  of  the  decrease  in  accuracy  when  observations  are 
not  available  in  the  immediate  vicinity  of  the  gun  emplacement, 

(b)  compare  the  difference  in  spatial  variability  derived  from  the  map 
winds  with  the  measures  of  spatial  variability  derived  directly  from  the  individ¬ 
ual  soundings, 

(c)  determine  if  topographic  effects  influence  the  areal  distribution  of 
measures  of  spatial  variability,  and 

(d)  determine  if  systematic  differences  exist  between  individual  single 
station  ballistic  winds,  and, where  such  differences  do  exist,  to  consider  the 
possibility  of  them  being  orogrsphically  induced. 


6.1  Areal  Distribution  of  the  rms  Differences  of  Analyzed  Ballistic  Winds 

Examples  of  fields  of  the  rms  difference  are  shown  for  Line-4  ballistic  winds 
in  Fig.  6-1  through  6-3,  for  Line  6  in  Fig.  6-4  through  6-6,  and  for  Line  10  in  Fig.  6-7 
through  6-9. 

As  would  be  expected,  the  magnitude  of  the  rms  differences  generally  increases 
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with  height.  The  rms  difference  for  the  v-component  of  the  ballistic  wind  exceeds  that 
for  the  u-component  for  all  levels  over  most  of  the  analysis  area. 

6.2  A  Comparison  of  Ballistic  Winds  Computed  from  Single  Station 

and  Analyzed  Winds 

An  important  characteristic  of  the  map  winds  is  illustrated  in  Fig.  6-10.  The 
rms  differences  of  the  two  components  of  the  ballistic  wind  based  on  the  CRAM  analyses 
for  Line  10  are  plotted  as  a  function  of  distance  from  the  trajectory  midpoint  along 
four  headings  as  indicated.  A  similar  statistic  was  computed  for  the  differences 
between  the  Ft.  Huachuca  single-station  ballistic  wind  and  ballistic  winds  derived 
from  five  stations  of  the  observing  network.  The  five  stations  compared  with  Ft. 
Huachuca  were  Benson,  Patagonia,  Parker  Canyon,  Bisbee,  and  Tombstone.  Their 
selection  was  based  on  all  available  1000  and  1400  MST  soundings.  The  rms  differ¬ 
ences  computed  from  observations  are,  in  general,  systematically  larger  than  those 
derived  from  the  analysis.  The  implication  here  is  that  if  several  observations  some¬ 
what  removed  from  the  firing  site  represent  the  only  ballistic  information  available, 
an  objective  analysis  of  these  observations  will  lead,  on  the  average,  to  an  estimate  of 
the  ballistic  wind  that  is  superior  to  the  ballistic  wind  derived  from  any  of  the  individual 
soundings. 

6.3  Systematic  Differences  Between  Single -station  Ballistic  Winds 

In  deriving  thu  rms  differences  for  station  pairs,  it  was  noted  that  there  were 
systematic  differences  between  the  ballistic  winds  derived  from  the  Benson  and  Tomb¬ 
stone  soundings.  In  particular,  for  Line- 10 ballistic  winds, for  those  cases  where  the 
wind  had  a  westerly  component  and  the  magnitude  of  the  wind  vector  exceeded  15  knots, 
the  u-component  of  the  Tombstone  ballistic  wind  exceeded  that  of  Benson  in  13  out  of 
15  cases.  The  average  difference  was  7.7  knots.  According  to  the  Student’s  T-test  for 
paired  comparisons  (Wadsworth  and  Bryan,  1960),  these  differences  are  significant  at 
the  1%  level. 

These  differences  may  be  traced  to  instrumentation  problems.  However,  an 
examination  of  the  trajectories  of  the  balloons  of  the  two  stations  provides  at  least  a 
clue  that  the  observed  differences  may  be  a  result  of  orographic  effects.  In  Fig.  6-11, 
Area  A  encloses  the  locations  of  the  Benson  balloon  at  the  time  it  was  passing  through 
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Zone  10  for  the  15  cases  studied.  Similarly,  Area  B  encloses  the  location  of  the  Tomb¬ 
stone  balloon  at  the  time  it  was  passing  through  Zone  10  for  each  of  the  15  cases. 

(There  are  two  such  areas  for  each  since  there  were  two  prevailing  wind  directions— 
northwest  and  southwest.)  In  most  instances,  the  Tombstone  balloon  at  the  important 
Zone- 10 levels  was  over  Sulpher  Spring  Valley,  in  contrast  to  the  associated  Benson 
balloons,  which  in  most  cases  were  well  to  the  west  over  higher  mountainous  terrain. 

There  was  also  noted  a  systematic  difference  between  the  Line- 10  single-station 
ballistic  winds  of  Ft.  Huachuca  and  Tombstone.  For  ballistic  winds  with  a  westerly 
component,  the  magnitude  of  Tombstones  winds  was  greater  than  that  of  Ft.  Huachuca 
in  10  out  of  11  cases.  The  average  difference  was  4.5  knots.  The  Student’s  T-test  for 
paired  comparison  indicates  the  differences  to  be  significant  at  the  5%  level. 

6.4  Some  Possible  Orographic  Effects  for  Low-level  Ballistic  Winds 

Detection  of  low-level  orographic  effects  was  rather  difficult  due  to  certain 
characteristics  of  objective  analysis  techniques.  In  general,  analysis  errors  are  a 
minimum  at  grid  points  near  the  observations  and  a  maximum  between  observations. 

To  reduce  the  analysis  error  maxima,  it  is  necessary  to  do  a  moderate  amount  of 
space  smoothing  (see  Section  3.3.5).  This  apparently  leads  to  a  reduction  in  the  esti¬ 
mates  of  the  rms  differences  at  grid  points  that  are  relatively  far  removed  from 
observing  sites.  Note, for  example,  the  secondary  maximum  in  the  rms-difference 
field  near  Nogales  in  Fig.  6-2.  Thus,  to  examine  the  rms-difference  fields  for  possible 
orographic  effects  in  the  lower-level  ballistic  wind  fields,  it  is  necessary  to  confine 
our  attention  to  rms  differences  at  grid  points  adjacent  to  observing  sites.  Figures  6-12 
and  6-13  show  the  rms  differences  for  Line-2  ballistic  wind  components  in  the  vicinity 
of  each  of  12  observing  stations.  A  minimum  in  the  rms  difference  for  the  v-component 
of  the  ballistic  wind  extends  up  the  San  Pedro  Valley,  which  includes  Benson,  while 
greater  variability  in  the  mountainous  terrain  can  be  noted  to  the  west,  south,  and  east 
of  Ft.  Huachuca.  The  rms  differences  for  the  u-component  of  the  Line -2  ballistic  winds 
yield  a  more  symmetrical  pattern,  with  only  a  slight  indication  of  a  minimum  in  the 
same  region  as  noted  above  for  the  v-component  differences. 
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6.5  Effect  of  Balloon  Drift  on  Single -station  Ballistic  Wind  Computations 

The  importance  of  accounting  for  the  drift  of  the  balloon  in  objective  analysis 
procedures  is  readily  apparent  from  an  examination  of  Figs. 6-7  through  6-9.  Super¬ 
imposed  over  the  analysis  of  the  field  of  rms  differences  in  this  figure  is  the  position 
of  the  Ft.  Huachuca  balloon  as  it  passed  through  Zone  10  for  those  cases  where  the 
wind  had  a  westerly  component  and  the  magnitude  of  the  wind  vector  exceeded  15  knots. 
It  can  be  seen  that  assuming  that  the  balloon  is  over  the  station  can  lead  to  significant 
errors  in  the  estimation  of  the  ballistic  wind. 
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Fig.  6-10.  Line-10  rms  space  differences,  u-  and  v-components. 


TABLE  7-2 

RESIDUAL  ERRORS  (in  Meters)  USING  INPUT  DATA  FROM 

Single  station  (ft.  huachuca)  and  cram  analyses 
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As  for  the  question  of  whether  the  single  station  observation  alone  results  in 
higher  residual  errors  than  those  given  by  using  values  obtained  from  a  CRAM  analysis 
with  one-half  the  data  as  input,  the  answer  is  not  readily  apparent  from  the  statistics 
for  the  five  cases  where  a  comparison  can  be  made.  The  average  absolute  and  rms 
residual  vector  error  is  slightly  low'er  when  the  Ft.  Huachuca  data  are  used  as  input. 

Examination  of  the  individual  component  error  reveals  that  the  range  errors  are 
lower  using  Ft.  Huachuca  data,  but  that  the  deflection  errors  are  lower  using  the  input 
from  the  CRAM  analysis  obtained  with  half  the  data.  Most  of  the  contribution  to  the 
lower  vector  errors  for  Ft.  Huachuca  data  is  due  to  the  lower  range  errors.  However, 
these  range  errors  have  been  empirically  corrected  (see  Section  4.2)  and  may  not  be 
as  reliable  as  the  deflection  error  values. 

Although  it  is  difficult  to  obtain  statistical  significance  with  the  small  number 
of  cases  in  this  sample,  the  difference  in  the  deflection  errors  data  was  subjected  to 
the  Student’s  T-testfor  paired  comparisons  (Wadsworth  and  Bryan,  1963).  The  result 
was  that  there  was  significance  at  the  15%  level.  The  value  for  the  level  of  significance 
is  the  same  as  the  probability  of  the  result  being  due  to  chance.  Thus,  there  is  a 
relatively  low  probability  that  the  better  result  given  by  the  CRAM  analyses  with  half 
the  data  is  due  to  chance. 

It  is  concluded  that  a  CRAM  analysis  which  uses  data  from  a  network  that  is  half 
the  density  of  that  designed  for  this  study  results  in  the  generation  of  significantly 
lower  residual  errors  for  the  deflection  component  of  the  error  than  those  provided 
by  data  from  a  single  station. 


8.0  PREDICTION  EXPERIMENTS 


To  develop  an  elaborate  system  to  forecast  ballistic  winds  in  mountainous 
terrain  is  an  enormous  undertaking.  In  this  study  the  objective  was  essentially  to 
assess  the  magnitude  of  the  problem  and  to  apply  a  few  very  modest  and  straight¬ 
forward  extrapolation  procedures. 

The  data  sample  used  for  these  experiments  consisted  of  CRAM  analyses  of 
ballistic  u-  and  v-wind  components  for  Lines  2,  4,  6,  8,  and  10  for  the  80  observation 
times  from  0600  MST  January  12  through  1400  MST  February  12,  1965.  Tne  observa¬ 
tions,  being  at  two-hourly  intervals  for  a  given  day,  allow  consideration  of  forecast 
lags  (" )  of  2,  4,  and  6  hr  with  valid  times  of  0800  (q=  2),  1000  (q=3),  1200  (q=4),  and 
1400  (q  =  5), 

There  were  three  principal  methods  used  to  forecast  ballistic  winds:  trend 
extrapolation,  linear  regression,  and  average  value.  Persistence  was  used  as  a  con¬ 
trol.  All  calculations  were  made  at  the  midpoint  of  the  firing  trajectory. 

(a)  Trend  extrapolation  There  were  two  types  of  trends  used;  a  2-hr 
trend  (2T)  that  calculated  the  difference  between  the  two  most  recent  consecu¬ 
tive  observations  at  the  trajectory  midpoint  and  extrapolated  that  difference 
forward  in  time,  and  a  4-hr  trend  (4T)  that  used  the  0600  and  1000  observations 
to  predict  for  1400.  The  following  equations  were  used: 

2T:  0(q)  =  0( q-r/2)  +  j  (0(q-r/2)  -  0(q-~^))  (8-1) 

for  r  =  2,  q  =  3,4,5 
r  =  4,  q  =  4,5 
r  =  6,  q  =  5 

4T:  0(5)  =  0(3)  +  (0(3)  -  0(1))  (8-2) 

where  0  is  the  extrapolated  value  of  the  observation  0  (either  u-  or  v-wi;(d 
component). 

(b)  Linear  regression  The  method  of  least  squares  is  used  to  fit  the 
0600,  0800,  and  1000  observations  to  extrapolate  for  1200  and  1400,  and  the 
0600,  0800,  1000,  and  1200  to  extrapolate  for  1400.  The  equation  is 

0(q)  =  a  +  bq  (8-3) 
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for  r  =  2,  q  =  4,5 
r  =  4,  q  =  5 

where  a  and  b  are  regression  coefficients  derived  from  q  =  1,  2,  3  or  q  =  1,  2,  3,  4.  ’ 

(c)  Averaging  The  average  value  of  the  past  and  present  observations  is 

computed  and  is  extrapolated  forward  as  a  constant  value. 

=  (8-4) 

8.1  Ballistic  Wind  Prediction 

Ballistic  wind  components  (u  and  v)  were  forecast  using  the  above  techniques. 

The  evaluation  measure  chosen  was  the  vector  rms  error.  The  results  are  shown  in 
Table  8-1  for  individual  lines  (2,  4,  6,  8,  10),  lags  (2,  4,  6  hr),  and  valid  times  (1000, 

1200,  and  1400).  Pooled  results  (all  hours)  are  shown  for  techniques  that  have  com¬ 
parable  data  samples.  Persistence  is  indicated  by  “P”  in  the  method  column.  The 
pooled  results  are  also  presented  in  Figs.  8-1  through  8-5.  It  can  be  seen  that  for  the 
lower  lines  (i.e.,  2  and  4)  the  averaging  technique  (AVE)  generally  yields  better  results 
than  the  2-hr  trend  extrapolation  (2T).  The  trend  extrapolation  errors  increase  with 
time  much  more  rapidly  than  do  the  AVE  errors  for  all  lines  (see  Figs.  8-1  through 
8-5),  especially  Line  2.  For  Lines  6;  8,  10,  the  errors  associated  with  2T  are 
smaller  than  AVE.  The  errors  associated  with  linear  regression  (LR)  were  roughly 
comparable  to  2T,  being  slightly  better  at  the  low  lines  and  slightly  worse  at  the  high 
lines . 

Of  interest  was  the  comparison  of  2T  and  4T,  in  which  2T  represents  the  fore¬ 
cast  for  1400  based  on  0800  and  1000  data  and  4T  represents  the  forecast  for  the  same 
time  based  on  f>600  and  1000.  For  the  lower  lines  (2,  4,  6),  the  4T  error  is  smaller 
than  the  2T  error,  but  for  Lines  8  and  10  the  reverse  is  true  although  the  differences 
are  small. 

Persistence  (P)  gives  the  smallest  errors  for  Lines  2  and  4  at  all  lags,  but  for 
the  higher  lines,  2T  is  generally  slightly  superior— the  exceptions  being  Lines  8  and  10 
for  t=6.  Figure  8-6  shows  a  comparison  between  P  and  2T  by  lag  and  line.  The  figure 
shows  that  2T  performs  best  at  Lines  4  and  6.  This  is  a  contrast  to  persistence,  which 
show's  an  increasing  rms  error  with  height. 
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TABLE  8-1 

BALLISTIC  WIND  PREDICTION  EXPERIMENT  RESULTS 
USING  CRAM  ANALYSES* 


(a)  Line  2 


Lag  (hrs) 

Method 

Valid  time 

1000  1200  1400 

All 

hours 

P 

2.7  3.5  3.3 

3.2 

o 

2T 

3.8  3.5  4.9 

4.1 

L 

LR 

3.9  4.3 

AVE 

3.0  4.5  4.2 

3.9 

P 

5.1  4.7 

4.9 

2T 

7.8  7.1 

7.5 

4 

4T 

6.3 

LR 

6.4 

AVE 

5.2  5.0 

5.1 

P 

5.4 

5.4 

6 

2T 

10.3 

10.3 

AVE 

5.5 

5.5 

(b)  Line  4 


Lag  (hrs) 

Method 

Valid  time 

1000  1200  1400 

All 

hours 

P 

2.9  3.8  3.1 

3.3 

2T 

2.9  4.1  5.2 

4.2 

2 

LR 

3.9  3.8 

AVE 

3.8  5.2  4.8 

4.6 

P 

5.4  4.7 

5.0 

2T 

5.9  6.7 

6.3 

4 

4T 

5.7 

LR 

5.6 

AVE 

6.1  5.9 

6.0 

P 

6.0 

6.0 

6 

2T 

7.6 

7.6 

AVE 

6.8 

6.8 

♦Values  are  root-mean- squart vector  errors  in  knots. 
Pooled  results  (all  hours)  ar  shown  for  techniques  which 
have  comparable  data  samples. 
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TABLE  8-1  (continued) 


(e)  Line  10 


Lag  (hrs) 

Method 

Valid  time 

'  All 
hours 

1000 

1200 

1400 

P 

5.0 

5.6 

4.1 

4.9 

O 

2T 

4.1 

4.3 

3.9 

4.1 

4 

LR 

5.3 

4.8 

AVE 

6.8 

9.1 

9.9 

8.6 

P 

9.8 

9.0 

mm 

2T 

9.5 

8.6 

mMm 

4 

4T 

8.9 

LR 

9.2 

AVE 

11.1 

12.1 

11.6 

P 

12.7 

12.7 

6 

2T 

13.8 

13.8 

AVE 

13.9 

13.9 

8.2  Artillery  Corrections  Based  on  Ballistic  Wind  Prediction  Experiments 

The  prediction  experiments  discussed  above  were  evaluated  on  the  basis  of  rms 
errors  where  the  verification  was  made  against  the  CRAM  analysis.  The  question 
may  be  raised,  quite  properly,  as  to  how  well  the  ballistic  wind  derived  from  CRAM 
represents  the  “true”  ballistic  wind.  This  can  be  answered  to  some  extent  by  con¬ 
ducting  an  evaluation  of  residual  errors  from  the  artillery  firings  based  on  forecast 
values  of  ballistic  winds,  temperatures,  and  density.  Because  trend  extrapolation 
yielded  results  as  good  as  any  of  the  other  methods  tried  (at  Line  10,  which  is  the 
only  line  for  which  the  artillery  firings  are  applicable),  it  was  chosen  for  evaluation, 
with  persistence  used  as  a  control.  Forecasts  were  based  on  both  CRAM  analyses 
and  the  Ft.  Huachuca  observation. 

Forecasts  were  valid  at  1000  or  1400  MST  on  days  when  artillery  firings  were 
conducted,  provided  that  the  meteorological  data  were  also  available. 

Two-hour  trend  forecasts  (2T)  were  made  using  0600  and  0800  data  for  a  valid 
time  of  1000,  and  1000  and  1200  data  for  a  valid  time  of  1400.  Four-hour  trend  fore¬ 
casts  were  made  two  ways  using  0800  and  1000  data  for  a  valid  time  of  1400  (2T,  r=  2), 
and  using  0600  and  1000  data  for  a  valid  time  of  1400  (4T). 
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Because  there  was  some  question  as  to  the  validity  of  the  range  corrections  (see 
Section  4),  the  evaluation  was  based  on  only  deflection  errors.  Unit  effects  were 
employed  with  the  forecast  parameters  in  the  same  manner  as  was  done  in  the  experi¬ 
ment  using  concurrent  data.  Table  8-2  lists  the  residual  deflection  errors  incurred  for 
various  lags  using  trend  extrapolations  and  persistence  based  on  CRAM  and  the  Ft. 
Huachuca  observation.  The  summarizations  in  Table  8-3,  based  on  a  limited  data 
sample,  indicate  the  following: 

(a)  Deflection  errors  based  on  persistence  are  smaller  for  CRAM  than 
the  Ft.  Huachuca  observations,  particularly  at  4  and  8  hr.  (Analysis  of  decay 
with  time  is  not  appropriate  here,  as  the  various  lags  involve  different  data 
subsets.) 

(b)  Trend  extrapolation  is  superior  to  both  CRAM  and  Ft.  Huachuca 
persistence  for  2-hr  lags,  but  not  for  4-hr  lags. 

(c)  A  4-hr  extrapolation  based  on  a  2-hr  trend  (2T,  r=4)  yields 
smaller  errors  than  a  straight  4-hr  trend  (4T),  but  not  smaller  than  persistence. 

The  time  decay  of  residual  errors  based  on  persistence  of  ballistic  winds  derived 
from  CRAM  analyses  was  examined.  Average  residual  deflection  errors  were  computed 
for  the  12  firings  conducted  at  1000  MST  for  0-,  2-,  and  4-hour  lags  and  for  the  12 
firings  conducted  at  1400  MST  for  0-,  2-,  4-,  6-,  and  8-hr  lags.  The  results  are  plotted 
in  Fig.  8-7.  Deviations  from  values  in  Table  8-3(a)  are  due  to  differences  in  sample 
size.  The  most  interesting  characteristics  of  Fig.  8-7  are  the  reductions  in  errors  for 
the  4-hr  lag  in  the  AM  firing  cases  and  for  the  8-hr  lag  in  the  PM  firing  cases.  In  the 
former,  the  error  decreases  from  130.0  meters  to  124.2  meters  when  going  from  2  to 
4  hr.  and  in  the  latter  the  error  decreases  from  146.9  meters  to  143.2  meters  when 
going  from  6  to  8  hr. 
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TABLE  8-3 

AVERAGE  RESIDUAL  DEFLECTION  ERRORS  (meters) 
FOR  PREDICTION  EXPERIMENTS 

(a)  CRAM  persistence  vs.  Ft.  Huachuca  persistence 


Lag  (r) 

All 

lags 

2 

4 

6 

8 

CRAM 

86.3 

81.5 

125.3 

109.3 

95.2 

Ft.  Huachuca 

88.9 

95.7 

127.3 

123.4 

103.4 

No.  of  cases 

17 

17 

9 

8 

51 

(b)  Trend  extrapolation  va.  persistence 


T 

=  2 

r 

=  4 

T 

=  4 

| 

Persist. 

Extrap. 

(4T) 

Persist. 

Extrap. 

(2T) 

Persist. 

CRAM 

77.4 

82.6 

144.5 

98.7 

116.0 

102.3 

Ft.  Huachuca 

87.5 

88.5 

169.9 

108.7 

135.4 

107.4 

No.  of  cases 

14 

14 

7 

7 

8 

8 

(c)  2T  (T  =  4)  vs.  4T 


2T 

4T 

Persist. 

CRAM 

140.3 

163.6 

113.7 

Ft.  Huachuca 

162.8 

196.3 

125.8 

No.  of  cases 

6 

6 

6 

Vector  rms  error,  knots 


14 


T,  hr 


Fig.  8-1.  Ballistic  wind  prediction  verification, 
Line  2. 
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Vector  rms  error,  knots 


Fig.  S-2.  Ballistic  wind  prediction  verification 
Line  4. 
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tor  rms  error,  knot 


Vector  rms  error,  knots 


Fig.  8-5 :  Ballistic  wind  prediction  verification, 
Line  10. 
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140 


afternoon  firings 


Fig.  8-7.  Average  residual  deflection  errors  based  on  CRAM 
persistence. 


9.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  analysis  technique  chosen  for  this  study,  CRAl*.  was  well  suited 
for  the  problem  at  hand  (i.e.,  development  of  techniques  for  correction  of  artillery 
firings),  and  produced  analyses  which  gave  a  reasonable  representation  of  the  wind, 
temperature,  and  density  fields  from  a  mesonetwork  of  twelve  rawinsonde  stations. 

The  analyses  that  were  generated  had  the  desirable  characteristics  expected  of  map 
winds,  i.e.,  relatively  smooth  contour  patterns  which  varied  slowly  with  time.  The 
application  of  the  method  to  three  dimensions  was  a  straightforward  extension  of 
existing  two  dimensional  CRAM  programs,  with  minor  modifications. 

Whe  .  applied  to  the  artillery  firing  situations  from  the  field  experiment,  meteoro¬ 
logical  parameters  derived  from  CRAM  yielded  smaller  residual  deflection  errors  on 
the  average  than  did  the  parameters  from  the  single-station  radiosonde  flight  (Ft. 
Huachuca).  The  average  deflection  error  for  CRAM  was  75.2  m  and  for  Ft.  Huachuca 
it  was  85.1  m,  based  on  19  firings.  The  average  range  error  for  CRAM  was  also  about 
1°  m  lower  than  Ft.  Huachuca,  although  the  absolute  magnitude  of  the  errors  was  sub¬ 
ject  to  question  due  to  a  possible  bias  in  the  firing  data. 

Time  and  space  variability  studies  were  conducted  to  study  the  effect  of  mountain¬ 
ous  terrain  on  the  variability  of  ballistic  winds.  The  time  variability  studies  indicated 
that  the  map  winds  derived  by  CRAM  gave  smaller  rms  time  differences  than  did  the 
single  sounding  winds.  The  rms  differences  at  Line  10  using  CRAM  were  as  much  as 
two  knots  less  than  the  Ft.  Huachuca  single  sounding  rms  differences.  The  same  was 
generally  true  for  the  lower  lines.  This  implies  less  decay  with  time  of  map  winds  than 
single  sounding  winds.  The  Line  10  rms  component  differences  at  the  trajectory  mid¬ 
point  increase  from  about  3—4  knots  for  2-hr  lags  to  10—12  knots  for  8-hr  lags.  This 

1/2 

is  more  rapid  than  the  a  =  2t  formula  would  predict.  This  variation  does  not  appear 

to  be  significantly  affected  by  change  of  sample  size  with  lag.  Differences  are  greater 

than  5  knots  for  a  4-hr  lag.  Line  4  field-average  rms  components  differences  increase 

1/2 

from  about  2  knots  at  2-hr  lag  to  about  5—6  knots  at  8-hr  lag  and  fit  c  =  2t  well. 

The  horizontal  range  of  time  variability  is  small  at  Line  10  but  is  larger  at  Line  2  and 
is  apparently  topographically  caused.  For  a  2-hr  lag  the  spatial  range  is  about  1/2 
knot  at  Line  10  and  about  one  knot  at  Line  2.  The  vertical  range  of  the  time  variability 
is  large.  There  is  a  general  increase  in  variability  with  height,  except  from  Line  2 
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to  4,  that  appears  to  be  caused  by  topography. 

The  space  variability  studies  suggest  that  in  a  situation  where  several  observa¬ 
tions  are  somewhat  removed  from  the  firing  site  and  they  represent  the  only  availaole 
source  of  ballistic  wind  information,  an  objective  analysis  (i.e.,  CRAM)  will  produce 
a  superior  ballistic  wind  estimate  than  any  of  the  observations  considered  individually. 

It  was  also  noted  that  there  was  a  significant  difference  in  the  u-component  of 
the  Line-10  ballistic  wind  between  Benson  and  Tombstone  when  the  winds  were  above 
15  knots,  with  the  Tombstone  u-component  averaging  7.7  knots  stronger  than  Benson’s. 
For  these  situations  it  was  found  from  the  balloon  trajectories  that  the  Bens  on-launched 
balloons  at  Zone  10  were  located  over  mountainous  terrain  while  the  Tombstone  balloons 
were  usually  over  relatively  flat  terrain  for  the  same  zone.  The  magnitude  of  the  Line- 
10  ballistic  wind  from  the  Tombstone  observation  was  also  4.5  knots  stronger  on  the 
average  than  the  Ft.  Huachuca  observation. 

Spatial  differences  at  low  levels  (i.e.,  Line  2)  were  also  attributed  to  terrain. 
Differences  in  the  v-component  toward  the  northeast  were  smaller  than  for  the  same 
distances  toward  the  more  mountainous  terrain  south  and  west  of  the  firing  site. 

In  determining  the  optimum  station  placement,  consideration  should  be  given, 
when  possible,  to  the  wind  climatology  of  an  area.  For  example,  more  observing  sites 
upstream  from  a  firing  range  and  fewer  downstream  would  allow  for  a  more  favorable 
spatial  distribution  of  balloons  when  they  reach  the  higher  zones,  e.g.,  Zone  10. 

Of  particular  importance  is  the  question  of  how  accurately  and  for  what  duration 
can  map  winds  be  forecast.  The  predictability  experiments  demonstrated  that  a  simple 
trend  extrapolation  at  or  near  the  trajectory  midpoint  will  not,  in  general,  produce 
ballistic  wind  superior  to  that  obtained  with  persistence.  However,  for  a  2-hr  lag, 
the  trend  extrapolation  does  give  smaller  residual  deflection  errors  than  persistence. 

For  14  cases  the  average  residual  deflection  error  was  77.4  m  for  a  2-nr  prediction, 
while  it  was  82.6  m  using  CRAM  persistence  and  88.5  m  using  Ft.  Huachuca  persistence. 
Th*3  comparison  between  CRAM  persistence  and  Ft.  Huachuca  persistence  showed  that 
the  CRAM  persistence  yielded  smaller  average  residual  deflection  errors  than  Ft. 
Huachuca  for  all  lags  tested  (2,  4,  6,  8  hr).  The  time  decay  characteristics  of  the 
deflection  errors  based  on  CRAM  persistence  showed  an  average  error  of  107.1  m 
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for  zero  lag,  120.4  m  for  a  2-hr  lag,  130.1  m  for  a  4-hr  lag,  146.9  m  for  a  6-hr  lag, 
and  143.2  m  for  an  8 -hr  lag. 

Many  of  the  conclusions  discussed  above  are  of  a  tentative  nature.  There  are 
several  areas  in  which  additional  investigation  should  be  carried  out  using  the  data 
that  have  already  been  edited,  error  checked,  and  put  on  magnetic  tape  for  input  to 
the  analysis  program.  For  example,  more  withheld  data  experiments  could  be  con¬ 
ducted  for  much  less  effort  relative  to  the  magnitude  of  the  data  editing  effort  that 
has  been  completed. 

Some  of  the  various  options  and  parameters  associated  with  the  CRAM  program 
should  be  tested  further.  Alternate  techniques  for  formulating  the  all-important  initial- 
guess  field,  such  as  use  of  a  prognosis, could  lead  to  significant  improvements  in 
determining  ballistic  winds. 

The  problem  of  predictability  needs  more  extensive  investigation.  A  simple 
approach  of  trend  extrapolation  was  only  successful  for  a  2-hr  lag.  The  prediction 
technique  described  in  this  report  made  no  distinction  among  the  various  zones  but 
rather  extrapolated  the  entire  line.  Additional  experiments  could  be  tried  in  which 
the  higher  zones  might  be  handled  differently  than  the  lower  ones.  The  results  of  the 
time  variability  studies  showed  that  the  variability  in  the  low  levels  was  much  smaller 
than  the  high  levels. 

An  alternate  method  for  predicting  the  higher  zone  winds  could  be  by  incorpora¬ 
tion  of  prognoses  from  existing  mid-tropospheric  NWF  models.  Still  another  possible 
approach  is  to  consider  both  space  and  time  in  extrapolating  the  wirds.  Consideration 
was  limited  to  only  time  extrapolation  in  the  present  study. 
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APPENDIX.  SYNOPTIC  SITUATIONS  FOR  DAYS  OF  BALLISTIC  WIND  ANALYSES 


Jan.  12,  1905 

At  sea  level,  a  large  nearly  stationary  high  located  over  eastern  Oregon  influ¬ 
enced  the  surface  winds  and  sky  cover  in  southeastern  Arizona.  At  0500  MST.one  hour 
prior  to  the  first  radiosonde  flight  for  the  small  network  of  stations  in  southeastern 
Arizona,  there  were  only  scattered  clouds,  light  southeasterly  winds  and  temperatures 
in  the  low  40’s  (*F).  Temperatures  warmed  to  the  low  60’s  by  1400  under  partly  cloudy 
skies  and  nearly  calm  surface  winds.  At  upper  levels,  a  weak  trough  was  oriented 
northeast-southwest  from  southern  Nevada  to  off  the  south  California  coast  at  0600, 
resulting  in  west  to  southwest  winds  aloft  with  speeds! about  30  knots  at  the  500-mb 

V 

level . 

Jan,  14,  1965 

The  sea-level  pressure  pattern  changed  little  from  that  of  Jan.  12.  The  high 
over  eastern  Oregon  moved  very  slowly  to  northern  Nevada  and  strengthened.  At 
0500  MST  the  surface  temperature  was  37°F  at  Ft.  Huachuca,  with  considerable  high 
cloudiness  and  no  wind.  The  temperature  climbed  to  the  middle  60’s  by  early  after¬ 
noon  and  conditions  remained  nearly  calm  under  high  clouds  that  allowed  sunshine 
through  (some  occasional  very  light  winds  from  the  east-southeast  were  reported). 

In  the  lower  and  mid  troposphere,  a  ridge  line  was  located  on  a  northwest-southeast 
line  from  a  high  center  off  the  Oregon  coast  to  southeastern  Arizona.  Thus,  the  winds 
were  very  light  and  variable  from  the  surface  to  the  mid  troposphere  (20,000  ft)  and 
only  increased  slightly  at  the  30,000-ft  level. 

Jan.  16,  1965 

The  high  pressure  system  over  the  western  states  on  the  12th  and  14th  remained 
strong  and  moved  only  to  Idaho  on  the  16th.  Gradient  flow  in  the  low  levels  over  south¬ 
east  Arizona  was  generally  from  the  east,  but  a  north-northeast  wind  at  10  knots  was 
reported  at  Ft.  Huachuca  at  0500  MST  with  a  temperature  of  45°F  under  cloudy  skies 
(high  clouds).  Afternoon  temperatures  rose  to  the  mid  60’s  near  Ft.  Huachuca,  but 
farther  west,  near  Tuscon,  temperatures  were  around  75°F. 

At  upper  levels,  a  strong  high  was  centered  almost  directly  over  the  sea -level 
high  in  Idaho  with  a  ridge  line  southeastward  over  southeastern  Arizona  causing  light 
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and  variable  winds  at  all  levels  to  30,000  ft. 


Jan.  18,  1965 

The  western  U.S.  continued  to  be  dominated  by  a  large  high -pressuresystem 
centered  in  Idaho,  resulting  in  generally  fair  skies  over  south-eastern  Arizona  all  day. 
Winds  at  the  surface  remained  nearly  calm  at  Ft.  Huachuca  for  both  the  1000  and  1400 
firings,  while  temperatures  climbed  rapidly  from  the  mid  30’s  at  0600  and  0800  to  the 
mid  70’s  by  the  1400  firing. 

At  upper  levels,  a  ridge  line  was  noted  from  western  Montana  southward  to 
central  Arizona,  this  ridge  caused  relatively  light  north  to  northwesterly  winds  over 
the  firing  area  during  the  early  morning.  The  winds  gradually  backed  to  more  westerly 
as  a  weak  trough  moved  eastward  across  the  west  central  states  during  the  day. 

Jan.  20,  1965 

The  high  pressure  system  centered  over  Idaho  for  5  days  was  still  there,  but 
had  weakened  considerably  and  a  trough  formed  from  central  Arizona  northwestward 
Some  showers  and  thunderstorms  occurred  along  the  trough  line  but  none  was  observed 
at  Ft.  Huachuca  through  the  day  to  the  1400  firing  time.  Temperatures  ranged  from 
the  upper  40’s  to  around  60°F  under  mostly  cloudy  skies,  and  surface  winds  were  south- 
southwest  around  10  mph  after  a  period  of  calm  from  about  0500  to  0800. 

At  upper  levels  a  weak  low  was  moving  southeastward  along  the  California/ 
Arizona  border  during  the  day,  causing  southwesterly  winds  about  20  knots  at  18,000  ft 
by  1400  MST. 

Jan,  22,  1965 

Southeastern  Arizona  was  under  the  influence  of  a  high-pressure  system  located 
over  Idaho  and  Nevada.  Once  again  surface  winds  were  very  light  or  calm  and  skies 
were  only  partly  cloudy.  Temperatures  in  the  low  30’s  at  0600  rose  to  the  mid  50’s 
by  the  1400  firing  time. 

At  levels  above  the  surface  a  deepening  trough  located  over  eastern  New  Mexico 
and  western  Texas  caused  northwesterly  winds  over  southeastern  Arizona  that  increased 
from  20—30  knots  at  about  10,000  ft  above  sea  level  to  50—60  knots  at  18,000  ft  and 
probably  higher  at  30,000  ft. 
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At  850  mb  (just  above  the  surface  of  the  firing  site)  the  winds  were  west-northwest 
at  5—15  knots. 

Jan.  25,  1965 

Southeastern  Arizona  was  dominated  by  the  circulation  around  an  intense  develop¬ 
ing  low  in  the  southern  plains  states  and  a  large  high-pressure  system  off  the  California 
coast.  Surface  winds  at  Ft.  Huachuca  which  were  light  at  0500  MST  increased  by  0800 
and  the  sky  was  obscured  by  blowing  sand  at  that  time.  Temperatures  ranged  from  an 
early  morning  low  around  35  to  a  high  of  about  50  at  the  1400  firing  time.  Skies  were 
mostly  clear  after  the  0800  sand  storm. 

At  upper  levels,  a  strong  trough  was  located  east  of  the  region  and  was  moving 
eastward  at  about  25  knots.  This  trough  caused  strong  northwesterly  winds  above  the 
Ft.  Huachuca  area  that  wfere  around  100  knots  at  the  30,000-ft  level  and  nearly  100  knots 
as  low  as  20,000  ft  at  0500  MST  . 

Jan.  27,  1965 

A  strong  high-pressure  system  moved  in  from  the  Pacific  Ocean  and  was  centered 
over  the  eastern  Arizona/Utah  border  during  the  day,  dominating  the  weather  over 
southeastern  Arizona  with  clear  skies  and  calm,  or  occasionally  light  easterly  surface 
winds.  Temperatures  ranged  from  around  25°F  at  0800  MST  to  the  low  40’s  by  the 
1400  firing  time 

A  trough  was  located  from  Central  Oklahoma  to  northern  Mexico  at  upper  levels 
resulting  in  north-northwesterly  winds  that  increased  with  elevation  to  around  40  knots 
at  18,000  ft  and  probably  stronger  at  25—30  thousand  ft  at  0500  MST.  The  winds 
diminished  somewhat  toward  mid-afternoon  as  the  trough  moved  farther  from  the  area. 

Jan.  29,  1965 

A  ridge  of  high  pressure  extended  southeastward  from  a  high  center  in  western 
Nevada  to  the  Ft.  Huachuca  area.  The  influence  of  ridge  conditions  resulted  in  clear 
to  partly  cloudy  skies,  calm  wind  conditions,  and  a  wide  temperature  range  from  24°F 
at  0800  MST  to  65°F  at  1400  MST. 

At  upper  levels,  southeastern  Arizona  was  in  a  "col”  area  with  generally  light 
winds  mostly  from  a  northerly  direction. 
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Jan.  31,  1965 

A  very  weak  high  system  was  centered  over  Arizona  causing  partly  cloudy  to 
clear  sky  conditions,  little  or  no  wind  at  the  surface, and  a  temperature  range  from  the 
low  30’s  at  the  time  of  the  0600  and  0800  radiosonde  flights  to  near  70“F  by  1400  MST. 

At  upper  levels,  westerly  winds  increased  with  increasing  height  and  were  grad¬ 
ually  shifting  to  more  northwesterly  through  the  day  and  strengthening.  Speeds  were 
about  25  to  35  knots  at  0500  MST  at  about  18,00  ft  and  had  increased  to  40  knots  at  the 
level  by  1700  MST. 

Feb.  2,  1965 

Throughout  the  day  a  stationary  front  lay  in  an  east-west  line  across  Northern 
Arizona  with  only  a  weak  pressure  gradient  over  southeastern  Arizona.  Thus,  winds 
were  light  and  variable  at  the  surface  and  skies  were  generally  cloudy  due  to  the 
proximity  of  the  front.  The  temperature  range  was  rather  small  from  around  50°F  at 
0500  to  the  low  60’s  by  the  1400  firing  time. 

At  upper  levels  a  weak  low'  was  centered  southwest  of  Arizona  just  off  the  south 
Pacific  coast  of  North  America  and  a  ridge  was  located  over  California.  Winds  were 
quite  light  at  levels  up  to  about  20,000  ft. 

Feb.  4,  1965 

Ridge  conditions  at  both  the  surface  and  aloft  prevailed  over  Arizona  all  day, 
with  winds  generally  calm  at  the  surface  and  light  northeasterly  aloft.  Temperatures 
ranged  from  near  30°F  at  0800  MST  to  the  mid  60’s  by  1400  MST. 

Feb.  6,  1965 

During  the  day,  a  developing  low-pressure  system  v/as  moving  southward  through 
eastern  Nevada  with  an  active  cold  front  extending  southward  from  the  low  center 
through  southeastern  California  and  off  the  Pacific  coast.  Showers  along  the  front  were 
advancing  eastward  across  Arizona,  but  had  not  yet  reached  the  Ft.  Huachuca  area  by 
the  1400  firing  time.  Skies  were  mostly  cloudy  through  the  day  near  Ft.  Huachuca  and 
surface  winds  were  increasing  from  the  south.  The  temperature  ranged  from  the  upper 
40’s  at  0600  MST  to  the  upper  50’s  in  the  early  afternoon. 

At  0500  MST  an  upper -level trough  extended  southward  from  central  Oregon  to  off 
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the  central  California  coast.  This  trough  progressed  slowly  eastward  during  the  day 
and  caused  upper-level  winds  over  southern  Arizona  to  be  from  the  west-southwest, 
and  later  to  be  from  the  southwest  with  increasing  speed.  Speeds  at  about  18,000  ft 
ranged  from  40  knots  at  0500  MST  to  about  60  knots  during  the  afternoon.  Veiy  likely 
the  winds  increased  with  elevation  above  18,000  ft. 

Feb.  8,  1965 

A  high  centered  over  Nevada  dominated  the  weather  over  the  southwestern  U.S. 

Winds  at  the  surface  were  light  and  variable  or  calm  all  day  and  skies  were  partly 
cloudy.  Temperatures  ranged  from  a  low  around  freezing  at  0500  MST  to  near  50°F  in 
the  early  afternoon. 

At  upper  levels  a  closed  low  was  centered  just  north  of  Tuscon  at  0500  MST, 
moving  slowly  eastward  during  the  day.  Winds  were  probably  from  the  southwest 
ahead  of  the  low  during  the  morning,  but  were  in  the  process  of  shifting  during  the  day 
as  the  low  passed  just  to  the  north  of  Ft.  Huachuca.  Speeds  ranged  from  20—35  knots  • 

at  the  18,000-ft  level. 

Feb.  10,  1965 

A  large  sprawling  trough  of  low  pressure  with  several  centers  was  located  to 
the  north  and  east  of  the  Ft.  Huachuca  area  over  Utah  and  New  Mexico  during  most  of 
the  day.  In  the  morning  hours,  0600  to  1200,  the  surface  winds  were  northerly  at 
5—10  knots,  but  backed  to  west-southwest  in  the  afternoon  and  increased  to  15—20 
knots.  Temperature  range  during  the  day  was  from  28°F  at  0500  MST  to  the  low  40’s 
during  the  afternoon. 

An  intense  upper -level  lowhad  moved  south-southeastward  into  northern  Arizona 
by  0500,  but  it  then  turned  east-northeastward  during  the  next  12  hours.  As  a  result, 
winds  were  quite  strong  (80—90  knots  at  about  18,000  ft)  throughout  the  day,  generally 
from  a  westerly  direction.  Wind  speeds  undoubtedly  increased  with  elevation  and  were 
probably  around  100  knots  near  the  30,000-ft  level. 

Feb.  12,  1965 

A  large  ridge  of  high  pressure  dominated  the  entire  western  quarter  of  the  U.S.  . 

(including  Arizona)  through  the  day.  Skies  were  nearly  clear  over  the  Ft.  Huachuca 
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area  and  calm  conditions  at  0600  and  0800  gave  way  to  west-northwest  surface  winds 
of  5—10  knots  for  the  1000  and  1400  firings.  Temperatures  ranged  from  a  low  of 
15°F  at  0600  and  0800  to  the  low  40’s  by  1400. 

At  ievels  above  the  surface,  a  rather  strong  ridge  of  high  pressure  that  extended 
in  an  arc  from  western  Canada  to  off  the  central  California  coast  resulted  in  north¬ 
westerly  flow  over  the  Ft.  Huachuca  area.  Speeds  increased  with  elevation  to  35—45 
knots  at  18,000  ft  and  probably  continued  increasing  to  the  30,000-ft  level. 
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